Quantum critical behavior in itinerant electron systems — Eliashberg theory and 
instability of a ferromagnetic quantum-critical point. 
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We consider the problem of fermions interacting with gapless long-wavelength collective bosonic 
modes. The theory describes, among other cases, a ferromagnetic quantum-critical point (QCP) and 
a QCP towards nematic ordering. We construct a controllable expansion at the QCP in two steps: 
we first create a new, non Fermi-liquid "zero-order" Eliashberg-type theory, and then demonstrate 
that the residual interaction effects are small. We prove that this approach is justified under two 
conditions: the interaction should be smaller than the fermionic bandwidth, and either the band 
mass tub should be much smaller than m — pf/vf, or the number of fermionic flavors N should 
be large. For an SU(2) symmetric ferromagnetic QCP, we find that the Eliashberg theory itself 
includes a set of singular renormalizations which can be understood as a consequence of an effective 
long-range dynamic interaction between quasi-particles, generated by the Landau damping term. 
These singular renormalizations give rise to a negative non-analytic q 3 ^ 2 correction to the static spin 
susceptibility, and destroy a ferromagnetic QCP. We demonstrate that this effect can be understood 
in the framework of the <f> A theory of quantum-criticality. We also show that the non-analytic q 3 ^ 2 
correction to the bosonic propagator is specific to the SU(2) symmetric case. For systems with 
a scalar order parameter, the q 3 ^ 2 contributions from individual diagrams cancel out in the full 
expression of the susceptibility, and the QCP remains stable. 

PACS numbers: 



I. INTRODUCTION 

Quantum-critical behavior in two-dimensional (2D) 
systems with continuous symmetry continues to attract 
substantial interest from the condensed-matter commu- 
nity. Near criticality, bosonic collective modes in ei- 
ther the spin or the charge channel (depending on the 
problem) are soft, and mutual feedback effects between 
bosonic and fermionic degrees of freedom lead to a rather 
peculiar behavior of both the fermionic and bosonic prop- 
agators. In this paper, we study in detail this behavior 
for two-dimensional (2D) systems. 

For systems with continuous symmetry, the dynamics 
of low-energy bosons is dominated by Landau damping, 
and the collective mode propagator dressed by a particle- 
hole bubble behaves as 
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(1.1) 



At £ = oo, the bosonic propagator becomes massless, 
signaling an instability towards a particular ordering. 
The dynamical exponent z, which measures how the fre- 
quency scales with momentum at criticality (u ~ q z ), is 
z = 3. 

Physically, the complexity of the problem resides in 
the presence of gapless fermions at the QCP. In ordinary 
QCP in localized electron systems, one deals with only 
one type of massless modes, namely a bosonic mode as- 
sociated with the fluctuations of the order parameter. In 
itinerant electron QCP, massless bosonic modes interact 



with conduction electrons, which are gapless at the Fermi 
surface, and this interaction affects both, electrons and 
bosons. One can still reduce the problem to interacting 
bosons by formally integrating the fermions out of the 
partition function. It was originally conjectured that 
this leads to a conventional </> 4 field theory with the bare 
propagator given by (|1.1|) . i.e., to a (j) 4 theory in an ef- 
fective dimension d + z. Since z = 3, the bosonic sector 
is above its upper critical dimension for d > 1, and the 
critical exponents have mean- field values. 

However, this description turns out to be oversimpli- 
fied by two reasons. First, it does not address the issue 
of what happens to the fermions at the QCP. It turns out 
that for fermions, the upper critical dimension is d^ = 3, 
such that in d = 2 the fermionic self-energy is singular 
and critically affects the behavior of low-energy fermions. 
Second, the 4> 4 theory for bosons is actually rather pecu- 
liar as the prefactors for the 4 and higher-other terms 
are determined by low-energy fermions, and are sensitive 
to to the behavior of these fermions near the QCP. 

These two arguments imply that at the QCP in itin- 
erant electron systems, fermions and bosons should be 
considered self-consistently and on equal footing. The 
key theoretical challenge in this context is to develop a 
controlled computational scheme to describe the correct 
behavior of electrons interacting with gapless bosonic col- 
lective degrees of freedom. 

The problem of fermions interacting with bosons with 
the propagator l|l.l|) at £ = oo was previously analyzed 
in the context of 2D fermions interacting with a singu- 
lar gauge field QEIQ, and was also applied to a gauge 
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theory of high T c superconductors |5| , and compressible 
Quantum Hall Effect @- Later, this problem was stud- 
ied in the context of 2D fermions near a ferromagnetic 
instability 0, H U EH , and, very recently, in the context 
of fermions near an instability towards a nematic-type 
ordering with angular momentum I = 2 [Til IT2I ITjjl uM . 
That last transition was argued by some studies [l5| to 
be relevant to the cuprates. Similar problems have been 
studied in the context of finite momentum spin |l6j and 
charge |17| ordering transitions in itinerant fermionic sys- 
tems. 

An analytic treatment of the problem was originally 
carried out by Altshuler, Ioffe and Millis (AIM) |j| 
for the interaction with the gauge-field. They showed 
that the fermionic self-energy scales as E(w m ) = 
*(|w m |) 2 / 3 cjQ^ 3 sign (u> m ) to second-order in perturbation, 
pointing to a breakdown of the Fermi liquid behavior. 
They estimated higher-order terms and argued that the 
w 2 /3 f orm f t ne self-energy survives to all orders in per- 
turbation. On the other hand, non-perturbative eikonal 
expansion and closely related approach based on 2D 
bosonization [3 yielded a different behavior, in 
which the fermionic Green's function decays exponen- 
tially with coordinate and time. This would be consis- 
tent with a divergent perturbative expansion for the self- 
energy. Altshuler and collaborators argued that this last 
result only survives in the artificial limit of a vanishing 
number of fermionic flavors N — > 0: at any finite N (in- 
cluding the physical case N = 1), the finite curvature of 
the fermionic dispersion prevents the perturbation series 
for the self-energy to become singular (see also Ref.polp. 

The discussion on the interplay between perturbative 
calculations and 2D bosonization re-emerged recently in 
the context of the quantum critical point for a Pomer- 
anchuk instability towards nematic ordering. Metzner 
and collaborators [13j and, very recently, Khveshchenko 
and one of us |2l| argued that E(w) oc ui 2 ^ 3 is the correct 
result at criticality, while Lawler et al 0] argued, based 
on 2D bosonization, that non-perturbative effects change 
this behavior. Adding to the controversy, Kopietz j^] ar- 
gued that higher-order corrections to the self-energy hold 
in powers of lo 2 ^ 3 \ logw| n , where the geometrical series of 
logarithms gives rise to an extra power of frequency, such 
that E oc uj 2 / 3 ~ a . 

In this paper, we re-analyze the problem. We con- 
sider a ferromagnetic QCP (z = 3 with spin SU(2) sym- 
metry), and a QCP towards nematic ordering, towards 
Ising-type ferromagnetism, and a gauge-field problem - 
these three last problems are mathematically equivalent 
and correspond to z — 3 and U(l) symmetry of the or- 
der parameter. We construct a controllable expansion 
at the QCP by creating a new, non Fermi-liquid "zero- 
order" theory by solving the set of coupled equations 
for the fermionic and bosonic propagators, while neglect- 
ing the vertex corrections as well as the momentum de- 
pendence of the fermionic self-energy. This procedure 
is often called Eliashberg theory because of its resem- 
blance to the Eliashberg theory for the electron-phonon 



interaction |23j. We analyze the residual interaction ef- 
fects using the new zero-order propagator instead of free 
fermions. We confirm an earlier result of that the 
residual interaction does not change the functional be- 
havior of the self-energy. We go beyond previous works 
and perform a careful analysis of the structure of the 
infra-red divergences in the theory. We analyze in detail 
the vertex corrections at various momenta and frequen- 
cies, the interplay between the Migdal approximation and 
Ward identities, the role of the curvature of the fermionic 
dispersion and the interplay between a direct perturba- 
tion theory for free fermions and an effective perturba- 
tion theory in which one expands around the Eliashberg 
solution. We also obtain the leading correction to the 
fermionic density of states (DOS). 

A generic condition for the validity of the Eliashberg 
theory is that bosons should be slow modes compared 
to fermions (i.e., for a fixed frequency, the bosonic mo- 
mentum should be larger than the fermionic one). Then 
fermions, forced by the interaction to vibrate at frequen- 
cies near the bosonic pole, are far from their own reso- 
nance and thus have a small spectral weight, giving rise to 
only a small correction to electron-boson vertex (Migdal 
theorem). Typical bosonic momenta in Eq. HI. If) scale as 
uj 1 / 3 and are obviously slower than free fermions whose 
momenta scale as u. The situation becomes less clear 
once the fermionic self-energy is included. We show that 
the correction to the static fermion-boson vertex is deter- 
mined by frequencies at which the fermionic self-energy 
is of order of a bare u>, and the static vertex correction 
is small, for a fermion-boson coupling smaller than the 
fermionic bandwidth. This coincides with the generic 
condition for the validity of the low-energy description 
since otherwise the physics is not restricted to the vicin- 
ity of the Fermi surface. 

This generic condition, however, is not a sufficient one 
at the QCP of spatially isotropic systems - we show that 
there are corrections to the Eliashberg theory which come 
from the scattering process in which one component of 
the bosonic momentum is near the bosonic mass shell, 
while the other is near the fermionic mass shell. We 
find that such corrections are dangerous if the expansion 
around the Eliashberg solution holds in powers of terms 
of order one. The same expansion around free fermions 
yields terms which formally diverge as powers of w -1 / 3 if 
one neglects the curvature of the Fermi surface 0, . 
We show, in agreement with [2j, that in this situation, 
the way to construct a fully controllable perturbation 
expansion around the Eliashberg theory at the QCP is 
to either assume that the curvature of the Fermi surface is 
large, or extend the theory to a large number of fermionic 
flavors, N. In this case, the self-energy diagrams with 
vertex insertions are all small, and the theory is under 
control. 

We emphasize that this smallness does not imply that 
the theory is in the weak coupling limit - the Eliashberg 
self-energy (given by the one-loop diagram) doesn't "feel" 
the curvature and diverges as a; 2 / 3 leading to non-Fermi 
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liquid physics at the QCP. Another exception is the pair- 
ing vertex, which doesn't feel the curvature as well, and 
is of order one. 

In the second part of the paper we show that there 
exists a third singular scattering process in which 
both fermionic and bosonic momenta vibrate near the 
fermionic mass-shell. This third process is qualitatively 
different from the first two scattering processes in which 
at least one component of the bosonic momentum is near 
the bosonic mass shell. 

We show that this process (which by virtue of scatter- 
ing near the fermionic mass shell is within the Eliashberg 
theory) gives rise to a non-analytic momentum expansion 
of the static vertex. We show that for a ferromagnetic 
SU(2) symmetric QCP, this non-analyticity eventually 
gives rise to a non-analytic and negative correction to 
the static spin susceptibility. This correction exceeds the 
q 2 term in and makes a ferromagnetic QCP unsta- 
ble. 

The issue of whether the QCP is internally stable has 
been the subject of numerous discussions in the recent 
literature. This work was pioneered by Belitz, Kirk- 
patrick and Vojta [24j who found that in a generic 3D 
Fermi liquid far from the QCP, the static spin suscepti- 
bility Xs(q) has a negative non-analytic momentum de- 
pendence, leading to a minimum of Xj 1 (?) located at 
some incommensurate momentum, rather than at q = 0. 
The same result was later obtained for 2D systems [26| . 
If one were to formally extend the Fermi liquid results 
to the quantum critical region, one would obtain that 
the continuous QCP becomes unstable [3(j • It was a pri- 
ori unclear, however, whether this extension procedure is 
justified, since the Fermi liquid behavior does not seem 
to survive as one approaches the QCP. 

To address this issue we explicitly compute the static 
spin susceptibility at criticality, and show that it is neg- 
ative and non-analytic at the smallest momenta. This 
implies that a ferromagnetic QCP is indeed unstable, as 
the Fermi liquid analysis suggests. We argue that the 
non-analyticity in Xs(q) is associated with the Landau 
damping term, which gives rise to an effective long-range 
dynamic interaction between quasiparticles both away 
from, and at the QCP. The singular fermionic self-energy 
at criticality only modifies, in not a very essential way, 
the functional form of the non-analyticity compared to 
that in a Fermi liquid. 

We also discuss the emergence of the non-analytic 
term in the static spin susceptibility in the framework of 
the Hertz-Millis-Moriya (HMM) (f> A theory of quantum- 
criticality |fj. As we said above, this theory assumes 
that there exists a regular expansion of the effective ac- 
tion in powers of the order parameter field <fi. We show 
that this is actually not the case, and the pre-factor of 
the </> 4 term is non-analytic and depends on the ratio be- 
tween typical momenta and frequencies. We show that 
this non-analyticity feeds back as a non- analytic \q\ 3 ^ 2 
correction of the quadratic term in <p. We study how the 
non-analyticity in Xs(l) affects the fermionic self-energy 



and show that it gives rise to self-energy terms larger 
that cj 2 / 3 (beginning at the three-loop order). The se- 
ries of such terms eventually leads to a breakdown of 
the Eliashberg theory for the fermionic self-energy at the 
energy scale related to the typical momentum scale at 
which Xs(l) becomes negative. 

We show that the non-analytic corrections to the 
bosonic propagator and the divergent corrections to the 
fermionic self-energy are specific to the SU(2) spin sym- 
metric case. For a nematic instability, as well as for a 
ferromagnetic instability in systems with Ising symmetry 
the dangerous terms in the bosonic propagator and the 
fermionic self-energy cancel out, and the QCP is stable. 

The paper is organized as follows. In Sec. II we dis- 
cuss the model. In Sec. Ill we present a quick analysis 
of the self-energies, justifying all at once, the Eliashberg- 
like treatment, and the need to include the curvature 
of the fermionic dispersion. In Sec IV, we discuss the 
Eliashberg theory near quantum criticality. In Sec. V, 
we analyze in detail the conditions one has to impose in 
order for the Eliashberg theory to be valid. This includes 
the computation of all vertex corrections, as well as the 
momentum-dependent self-energy at the two-loop level. 
We also compare the results obtained by strict perturba- 
tive expansion using free fermions, and the Eliashberg- 
type calculations. 

In Sec. VI we address the issue of the stability of a fer- 
romagnetic QCP. We revisit the scaling arguments that 
a ferromagnetic QCP must be stable in dimension D > 1 
and show that the prefactor for the </> 4 term is actually 
a nonanalytic function of the ratio of frequency and mo- 
mentum. We show that this non-analyticity feeds back as 
a non-analytic correction to the static spin susceptibility. 
We explicitly compute the momentum-dependent term 
in Xs(q) at the two-loop order, both in the Fermi liquid 
regime away from a ferromagnetic QCP and at criticality. 
In both cases, we find that the dominant term at small q 
is negative and non-analytic. We also show that the in- 
stability of a ferromagnetic QCP can be also detected by 
computing the fermionic self-energy which at the three- 
loop order acquires extra singular terms because of the 
singularity in the static susceptibility. In Sec. VII, we 
evaluate the charge susceptibility at two-loop and three- 
loop orders and show that it remains analytic - all non- 
analytic contributions from individual diagrams cancel 
out. Finally, in Sec. VIII, we present our conclusions 
and discuss the consequences of the instability of a con- 
tinuous QCP towards ferromagnetic ordering. Technical 
details are presented in the Appendices A-F. 

A short version addressi ng p art of the results of section 
VI has been presented in |31j. 



II. THE MODEL 

The model we consider describes low-energy fermions 
interacting with Landau-overdamped collective bosonic 
excitations which are either gapless by symmetry reasons 
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(as it is the case for the interaction with a gauge field), 
or become gapless at the quantum critical point (for the 
nematic and the ferromagnetic problems). 

The underlying lattice models may be quite different 
for these three cases, but the low-energy models are very 
similar, the only difference being that in the case of the 
ferromagnetic QCP the gapless bosonic excitations are 
in the spin channel, whereas in the nematic case and the 
gauge field problem they are in the charge channel. 

The general strategy to derive the low-energy model is 
to start with a model with fermion-fermion interaction, 
assume that there is only one low-energy collective de- 
gree of freedom near the QCP, decouple the four-fermion 
interaction term using the critical bosonic field as an 
Hubbard-Stratonovich field, and integrate out of the par- 
tition function all high-energy degrees of freedom, with 
energies between the fermionic bandwidth W and some 
cutoff A [llll^. 

If this procedure was performed completely we would 
obtain a full Renormalization Group treatment of the 
problem. Unfortunately, there is no controllable scheme 
to perform this procedure. It is widely believed, though 
that although the propagators of the remaining low- 
energy modes possess some memory of the physics at high 
energies, the integration of high-energy fermions does not 
give rise to anomalous dimensions for the bare fermionic 
and bosonic propagators in the low-energy model. In 
practical terms, this assumption implies that the bare 
propagator of the relevant collective mode is an analytic 
function of momentum and frequency, and the fermionic 
propagator has the Fermi liquid form: 



(2.1) 



where Zq < 1 is a constant, and e k is the renormalized 
band dispersion. Near the Fermi surface, 



e k = v F k± 



2m B 



Here k is the momentum deviation from kp, the parallel 
and perpendicular components are with respect to the 
direction along the Fermi surface at kp, mj is the band 
mass, the Fermi velocity vf = kp/m, and for a circular 
Fermi surface one has m = tub- 

One can then re-cast the original model of fermion- 
fermion interaction into an effective low-energy fermion- 
boson model. Consider for dcfiniteness that the system 
is close to a ferromagnetic QCP. Then the low-energy 
degrees of freedom are fermions (with the propagator 
given by H2.1|l 1 and long- wavelength collective spin exci- 
tations whose propagator (the spin susceptibility) is an- 
alytic near q — and Q = 0: 



Q 



An 2 + o{q 4 ,n 4 )' 



(2.3) 



Here A is a constant, and £ is the correlation length, 
which becomes infinite at the QCP. We prove in the next 



section that the f2 2 term does not play any role in our 
analysis, and we therefore neglect it from now and ap- 
proximate the above bare propagator by the static one 
Xs,o(q)- The model can then be described by the phe- 
nomenological spin-fcrmion Hamiltonian: 



k,a q 

+■9 E cl a a al3 c k+q ^ ■ S q . 



(2.4) 



Here S q with q < A/vp are vector bosonic variables, and 
g is the effective fermion-boson interaction. For conve- 
nience, we incorporated the fermionic residue Zq into g. 

To illustrate how this effective Hamiltonian can, in 
principle, be derived from the microscopic model of inter- 
acting conduction electrons, we consider a model in which 
the electrons interact with a short range four-fermion in- 
teraction U(q) and assume that only the forward scatter- 
ing is relevant (U(0) — U): 



H 



k . a 



+ \zZ U zZ C fe,a C fc+9," C fc'/3 C fc'-9:/3; ( 2 - 5 ) 
q k,k',a,/3 

In this situation, the interaction is renormalized inde- 
pendently in the spin and in the charge channels |33| . 
Using the identity for the Pauli matrices a a g ■ cr~ ( s = 
+ 28 a s5f3 1 , one can demonstrate |3i| that in each 
of the channels, the RPA summation is exact, and the 
fully renormalized four-fermion interaction {q) is 

given by: 



Sa^Spe [ - + Gp 



1 



(2.2) 

where tr^ are Pauli matrices (a = x,y,z), and 



yp ~ 2 1 - un ( q ) ' 2 1 + un ( q ) ' 



(2.6) 



(2.7) 



with n(g) = -£(1 - a 2 (q/k F ) 2 ), a > 0. 

For positive values of U satisfying mU /2tt w 1, the in- 
teraction in the spin channel is much larger than the one 
in the charge channel. Neglecting then the interaction 
in the charge channel, we can simplify the Hamiltonian 
(1231) : 



H = Yl ekC k,a C k,* + |Eg U ef f(q) 



k,oc 

x E 



cl,a a af}Ck+q,f3 ' Cy^gCV -q,S- 



(2.8) 



where U eff (q) = (l/2)U 2 U(q)/{l + Utt(q)). Performing 
a Hubbard-Stratonovich decomposition in the three fields 
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k. 



c) 



FIG. 1: a) Polarization bubble b) One-loop fermionic self- 
energy c) Two-loop fermionic self-energy. 



S q , one recasts (|2.8|) into Eq. H2.4[l with: 



r a 

Xo 
9 

It 2 



U- 

2„ 



g ¥ Xo = (U/2)k% 



(2.9) 



( 



_2tt_ 

mil 



The QCP is reached when mU/2ir = 1, i.e., £~ 2 = 0. 
This coincides with the Stoner criterion for a ferromag- 
netic instability [34|. 

We emphasize that the bosonic propagator in Eq. (|2.3|) 
does not contain the Landau damping term. This is be- 
cause we only integrated out the high-energy fermions, 
whereas the Landau damping of a collective mode of en- 
ergy f2 comes from fermions of energy to < f2, and can 
only be generated within the low-energy theory. The dy- 
namics of both the bosonic fields S q and the fermionic 
c and is determined self-consistently by treating both 
fluctuations on equal footing. 

To put under control the computations carried out 
later in the paper, it is necessary to extend the model 
by introducing N identical fermion species, while keep- 
ing the SU(2) spin symmetry. The Hamiltonian (|2.4() 
can then be rewritten as: 

H s f = Hf + Hb + Hi n t where 

H b = X^o(<?) S 9 ' S -q 

q 

Hint = g C \^,a a al3Ck+q,j,(3 ' Sq, (2.10) 

k,q,j,a,/3 

where the index j = 1 . . . N labels the fermionic species. 

We use the spin-fermion Hamiltonian of Eq. (|2.10() 
as the starting point of our analysis. In the case of a 
QCP in the charge channel, or a ferromagnetic instability 
with Ising symmetry, the bosonic vector field S becomes 
a scalar field designated as </>. The interacting term is also 
modified, the Pauli matrices being replaced by 5 a p for the 
interaction with charge fluctuations, and by in the 
Ising case. The corresponding interaction Hamiltonians 
are: 



H 



Charge 

int 

Ising 



Et Z J. K^- 11 ) 

k,q,j,a,0 C k,j,a a a/3 Ck +<ld,l 3( Pq 



III. DIRECT PERTURBATION THEORY 

In this section we compute the fermionic and bosonic 
self-energies for the model presented in Eq. (|2.4I) using a 
perturbation expansion around non-interacting fermions. 
Our goal here is three-fold: to relate the Landau damping 
coefficient to the fermion-boson coupling constant g, to 
distinguish between E(w) and S(fc) and to demonstrate 
the importance of the curvature of the Fermi surface. 

We evaluate the self-energy in this section up to two- 
loop order. We verified that to this order, there is no 
qualitative difference between the quantum critical point 
in the spin or in the charge channel. We then restrict our 
presentation to the spin-fermion model near a ferromag- 
netic QCP. 



A. Bosonic self-energy: the Landau damping term 

The full bosonic propagator depends on the self-energy 
H(q, ri) according to: 



Xo 



^- 2 + q 2 +U(q,n) 



(3.1) 



At the lowest order in the spin-fermion interaction, the 
bosonic self-energy is given by the first diagram in Fig. 
^ and reads: 



II(q,fi) = 2Ng J G(k, W ) G(k+q, 



to+n). (3.2) 



The curvature of the fermionic dispersion does not af- 
fect much the result of this computation as it only leads 
to small corrections in q/rriBVF ■ Neglecting the quadratic 
term in the fermionic propagators, we introduce the an- 
gle 9 defined as £k+q = £k + VFq cos 9 and perform the 
integration over et, which gives us: 



n(q,fi) 



. Ngm 



oc 

271 



1 



9{u)) 



Nmg 



x / d9 

o til — VFqcoa9 

101 



TT ^{v F qY + ft 2 



(3.3) 



At the QCP, the bosonic mass-shell corresponds to the 
region of momentum and frequency space for which the 
terms in the inverse propagator are of the same order, i.e. 
near a mass shell q and fl satisfies II(q, ft) ~ q 2 . It follows 
that, at the QCP, near the bosonic mass shell, vpq/U ~ 
VFimgv 2 ? /fi 2 ) 1 / 3 ^ 1 at small enough frequency, so that 
VFq is the largest term in the denominator of H(q,il). 
The expression of the bosonic self-energy then reduces 
to: 



n(q,fi)=7 



101 



(3.4) 
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where 7 = 



Nmg 

TTVf 



This expression describes the Landau damping with a 
pre-factor depending on the fermion-boson coupling con- 
stant. This term is larger than a regular 0(fl 2 ) term, and 
fully determines the dynamics of the collective bosonic 
mode. 



and critical bosons leads to a breakdown of the Fermi- 
liquid behavior: at low energies, the w 2 / 3 term in i|3.7|) is 
larger than the bare u> in the fermionic propagator. This 
obviously makes one wonder if higher order insertions 
lead to even more singular contributions. 

We next compute the one-loop self-energy E(k), given 
by: 



B. One-loop fermionic self-energy 

We now turn to the fermionic self-energy, right at the 
QCP, where £ _1 = 0. To the lowest order in the interac- 
tion, the fermionic self-energy contains one bosonic line, 
as represented in Fig^j, and its analytic form writes: 



d 2 q dfl 
(2tt) 3 



Go(k + q,co + tt)x(q,fi). 



(3.5) 

The superscript "free" implies that we use the free 
fermionic Go(fc,c<j) in the integral for the self-energy. 

In a direct perturbation theory in g, one would have to 
use the bare form of the bosonic propagator, Eq. I|2.3[l . 
which leads to X(u) oc logw. However, this result is use- 
less, as we already know that the Landau damping term 
completely overshadows a regular frequency dependence 
in (|2.3[l . It makes more sense then to estimate the per- 
turbativc self-energy using the full bosonic propagator 
Eq. (|3.1(l . This is not a fully self-consistent procedure, 
but we use it here to estimate the functional forms of the 
self-energy at various orders in perturbation around free 
fermions. 

It is instructive to distinguish between £(k F ,c<j) = 
S(w) and £(k,0) = £(k). 

Substituting the renormalized bosonic propagator with 
the Landau damping term into Ij3.5|l . the frequency- 
dependent self-energy reads: 



,f ?o / dfl qdq d6 tttt 

(2w) 3 J q 2 +^f 

1 



i(u> + f2) — vpq cos ( 



(3.6) 



Here 8 is the angle between k F and q, and we linearized 
the fermionic dispersion. Evaluating the integral over 
the angle, and using that the typical internal bosonic 
momentum q ~ 5! 1 / 3 is much larger than fi ~ lj, we 
obtain 



~i free 



M = 



3.9 



where 



lu 



2tt 2 vf 



, , 1/3 < , 2 / 3 



3V3ff 3 

8lT 3 V F ^ 



dCl 



dqq 



3V3g 2 
8n 2 Nmv 2 F ' 



(3.7) 



(3.8) 



This result has been first obtained in Ref. It shows 
that in D = 2, the interaction between bare fermions 



E^ ree (k) = 



WfJ dn qdq dd J7^m. 



iQ — €k — VFq cos 9 



3.9 

(2ny 



dQ dq 



q \n\ 



f + 7 M (n 2 + (v F q) 2 ) 3 / 2 



.(3.9) 



One can make sure that the integral is infrared- 
convergent, i.e. S frcc (ek) oc ek, with an interaction de- 
pendent prefactor, which also depends on the upper cut- 
off of the theory, A. This suggests that the momentum 
dependent part of the fermionic self-energy is regular at 
the QCP and only leads to a finite mass renormalization. 



C. Two-loop fermionic self-energy 

We next calculate the contribution to the fermionic 
self-energy S(w) from diagrams at the two-loop level. For 
illustrative purposes, we consider the diagram presented 
in Fig. ^i, which writes: 

£2° C M~3 2 / diJ 1 d 2 q 1 [ dw 2 rf 2 9 2 G(k F + q 1 ,w + cj 1 ) 



x G(k F + qi + q2, w + toi + lu 2 ) G(k F + q 2 , w + uj 2 ) 
x X s (qi,wi) x s (q2,w 2 ), 

(3.10) 

where we use the full bosonic propagator, the free 
fermionic one, and we restrict ourselves to the frequency 
dependence. 

We first compute this integral expanding the disper- 
sion of the internal fermions to linear order, since the 
quadratic term was not significant in the computation of 
the one-loop bosonic and fermionic self-energy. Choosing 
the x axis along the external k = k F and integrating over 
qf and qf, we are left with: 



£ 2 ree M 



9 

V%LO 



duj\dq 



7l"i| 
lly 



du>2<lq: 



2y 



"1 lly 



7l"2| 

12 y 



) 2 / 3 ,„l/3 



(3.11) 



where too is defined in l|3.8JI . At low-energy, this two- 
loop self-energy diverges faster than the one-loop self- 
energy obtained in (|3.7|l . Estimating higher-order di- 
agrams, we find that they form a series in powers of 
(wq/w) 1 / 3 , such that the perturbative expansion around 
free fermions breaks up at u ~ loq. This result is in 
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line with the one obtained by Lawler et al. |l2| us- 
ing a two-dimensional bosonization scheme. The scale 
luq is related by luq — Vf/xq to the spatial scale at 
which the equal time fermionic Green's function G{x), 
obtained from bosonization, begins decaying exponen- 
tially (G(x) oc e-^l*^ 1 ' 3 ). 

However, the divergence of the perturbation theory can 
be cured once the curvature of the fermionic dispersion is 
included, as we now show. We re-evaluate the two-loop 
self-energy (I3.10|) . using now the full fermionic disper- 
sion, Eq. (|2.2(l . After integrating over the momentum 
component qf and gf , one has: 



dtoidq 



qf + 



7l"i| 
Hy 
1 



doj 2 dq: 



■Mi 



7 1^2 | 

I2y 



1 —1 

m g , 2 



(3.12) 



This result agrees with 0. We see that, when the cur- 
vature of the fermionic dispersion is included, that the 
two-loop self-energy turns out to be small compared to 
its one-loop counterpart, at low energy. In a separate 
study |2l|, one of us (A.C.) and D. Khveshchenko recon- 
sidered the bosonization procedure in the presence of the 
curvature and obtained the same results as in (|3.12(l . 



D. Summary 

As a conclusion, this first approach suggests that both 
the fermionic and the bosonic self-energies are important 
at the QCP. The bosonic self-energy sets the dynamics of 
the bosons, while the fermionic self-energy is non- analytic 
and parametrically larger than the bare u term at low 
energy, which implies a breakdown of the Fermi-liquid 
behavior at criticality. 

We also found that only the frequency-dependent part 
of the self-energy matters, the momentum-dependent one 
only leads to a regular renormalization of the effective 
mass. Finally, we found that the curvature of the Fermi 
surface plays an important role in regularizing the per- 
turbation expansion. 

The full account of these effects cannot be obtained 
from this simple analysis and one has to develop a con- 
trollable way to treat the bosonic and fermionic self- 
energies on equal footing. Since we found that only 
the frequency-dependent S(w) is relevant, a way to pro- 
ceed is to verify whether an Eliashberg-like theory, sim- 
ilar to the one developed in the context of phonon 
superconductivity |2^|, may be such a controllable ap- 
proximation. 



IV. ELIASHBERG THEORY 

The Eliashberg procedure allows to compute the 
fermionic self-energy £(w) and the bosonic polariza- 
tion II(<7, SI), by solving the self-consistent set of cou- 
pled Dyson's equations, neglecting all contributions com- 
ing from the vertex corrections and the momentum- 
dependent fermionic self-energy. 

Specifically the Eliashberg theory follows three steps: 

• neglect both the vertex corrections and the momen- 
tum dependent part of the fermionic self-energy, 
i.e., approximate: 



S(k,w n ) = S(w„) 
5Tot = 9 + Ag = , 



(4.1) 



• construct the set of self-consistent Dyson's equa- 
tions: 



G (k, tj n ) = i(J n — vp(k — ftp) + iT,(uj n ) 
X s (g,Sl m ) — 



Xo 



r 2 +q 2 + U(q,Q m ) 



(4.2) 



with the following fermionic and bosonic self- 
energies: 



«E(w„) 



k, uj„ 



-1 / \ <l> v mf~~\ 

Xo n(g,S2 m ) = vyi jw 



(4.3) 



The fermionic Green's functions in l|4.3|) are 
full (they are represented diagrammatically by a 
straight line) and Xs{q, ^m) is the full bosonic prop- 
agator (represented by a wavy line). 

• check a posteriori that the neglected terms Ag and 
E(k), are all parametrically small. 

The evaluation of the momentum integral for the 
fermionic self-energy in the Eliashberg theory requires 
care. Since fermions are faster than bosons, the leading 
contribution to S(w) is obtained if one integrates over the 
momentum component transverse to the Fermi surface 
only in the fermionic propagator and sets this component 
to zero in the bosonic propagator (this implies that the 
momentum integral is factorized) . One can show that the 
corrections that arise from keeping the transverse com- 
ponent of momentum in the bosonic propagator are small 
to the same parameter as Ag/g and should therefore be 
neglected, as keeping them would be beyond the accu- 
racy of the theory. The factorization of the momentum 
integration is what distinguishes the Eliashberg theory 
from the FLEX approximation. In the FLEX approx- 
imation, one also neglects vertex corrections, but does 
not factorize the momentum integral in (|4.3() . 

The evaluation of the bosonic and fermionic self- 
energies within the Eliashberg theory is presented in Ap- 
pendices El an d DH We list here the main results. 
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At large but finite correlation length £ and for a 
bosonic momentum and frequency satisfying vpq 3> 
E(fi) we obtain: 



U(q,Q) =T~Y ,and 
E(w) = AcjF (j^ 3 ) , 



(4.4) 



where F(ac < 1) = 1 + O(x), and F(x > 1) = ^2;~ 1/3 - 
The parameter 7 is the same as for free fermions, 
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Nmg 
TTVf ' 



and A is the dimensionless coupling 
35 



A 



ff 2 Xo 



(4.5) 



(4.6) 



At finite £ 1 and vanishing u>, the self-energy has a Fermi 
liquid form: 



E(u;) = Xlo. 



(4.7) 



The Fermi liquid theory is stable, and the low-energy 
quasi-particles have a finite effective mass: 



= m(l + A). 



(4.8) 



The effective mass diverges proportionally to £ in the 
vicinity of the QCP. 

At w > (7£ 3 ) -1 , however, the system is in the 
quantum-critical regime. Here the Fermi liquid theory 
breaks down in the sense that the quasi-particles damp- 
ing becomes comparable to its energy. We have: 



E(w) = w 1/3 M 2/3 signH, 



(4.9) 



where uj = 3\/3g 2 / (8ir 2 Nmv 2 F ) is the same as in l|H.8|l . 

At the QCP, = 0, the region of Fermi- liquid behav- 
ior collapses, and the lo 2 / 3, dependence of the self-energy 
extends down to the lowest frequencies. The expression 
for E(w) is valid for all frequencies below the cutoff A. 
However, only frequencies u> < loq are actually relevant as 
at higher frequencies the system behaves as a nearly ideal 
Fermi gas. Note that the curvature of the fermionic dis- 
persion is unimportant here and only accounts for small 
corrections containing higher powers of frequencies. 

Comparing l|4.4|) and l|4.9|l with (|3.4|) and (|3.7|l . we 
see that the self-energies in the Eliashberg theory coin- 
cide with the one-loop perturbative results around free 
fermions. This arises from the fact that the momentum 
integration is factorized in the Eliashberg theory, and the 
full fermionic propagator appears in both self-energies 
only via the fermionic density of states (DOS): 



N(w) 



dei 



' LU + E(u>) — 6fe 



(4.10) 



This DOS reduces to N(lu) = sign{uj), independently on 
the self-energy: it remains the same as for free fermions. 



Note that Eq. (|4.4() for the bosonic self-energy is only 
valid as long as the interplay between the external mo- 
mentum and frequency is such that vfQ^ E(f2). In the 
opposite limit, the vertex corrections cannot be neglected 
as we argue in the next section. 

For further analysis of the corrections to the Eliash- 
berg theory, it is then convenient to rescale N out of the 
formulas for E(w) and for H(q, f2). This can be done by 
rescaling m and fcj? leaving vf intact: 

m -> m/N k F -> k F /N. (4.11) 

We emphasize that the fermionic self-energy E oc cj 2 / 3 
and the Landau damping term in the bosonic propagator 
do not contain N after rescaling, and therefore must be 
included into the new "zero-order" theory about which 
we then expand using X/N as a small parameter. This 
zero-order theory includes: 

Xo 




1 



•7|fi|/9 



(4.12) 



i(u> + E(w)) - e k ' 



At the QCP, this theory becomes completely parameter- 
free once we measure the frequency in units of uiq and 
rescale the bosonic momentum. 

We show in the next two sessions that the corrections 
to the Eliashberg theory are small in powers of two pa- 
rameters, one of which vanishes in the limit N — > 00, 
the other being parametrically small. Accordingly, the 
Eliashberg theory becomes exact at N — 00, however the 
last limit has to be taken at a finite curvature. 

For the SU(2) ferromagnetic case, we will show below 
that Eq. (|4.12l) is incomplete. There are extra contri- 
butions to both xility an d E(a>) which also belong to 
the Eliashberg theory, showing up at three-loop level and 
higher for the fermionic self-energy. These terms, how- 
ever, appear due to rather specific reasons related to the 
presence of an effective dynamic long-range interaction in 
the theory and all cancel out for the gauge field problem, 
at a nematic QCP and for an Ising ferromagnetic QCP. 
We consider these extra terms in the next section and 
proceed here without taking them into consideration. 

We can reformulate the Eliashberg theory by intro- 
ducing the following effective Lagrangian describing the 
fermion-boson interaction: 

L = L F + L B + L int , with 

Lb =T^2S q , n X~ 1 (q,£ln)S-q-n, 

q,n 

L F =T ^ C fej\n,a G ' _1 ( fc ' Li; n) C fcj 1 n,"' 
k 7 n,a 7 j 

Lint =gT 2 S gj7n C k j na (T a 0Ck-q,j,n-m,l3) (4-13) 

where n, m number Matsubara frequencies, a, (3 are spin 
indices, and j is a flavor index. The upper limit of the 
frequency summation is the cutoff A. 
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We emphasize that there is no double counting in 
the bosonic propagator l|4.13|) . The integration of high- 
energy fermions, above the cutoff A leads to the mo- 
mentum dependence of the static bosonic propagator, 
whereas the interaction at frequencies below A gives rise 
to the Landau damping, without affecting the static part. 

To summarize, the Eliashberg-type theory at the QCP 
contains a non-analytic fermionic self-energy which scales 
as a; 2 / 3 and breaks down the Fermi liquid description 
of fermions. At the same time, the bosonic propagator 
is regular - the only effect of the interaction with low- 
energy fermions is the appearance of a Landau damping. 

We will see below that for a charge QCP ( and for a 
spin QCP with a scalar order parameter) this Eliashbcrg 



theory is entirely stable. For a ferromagnetic SU(2) sym- 
metric QCP, the Eliashberg theory has to be extended, as 
we will see in Sec VI, to incorporate extra singular terms 
associated with the existence of long-range dynamic in- 
teraction between quasi-particles. 

We emphasize that even for a charge QCP the fully 
rcnormalizcd bosonic susceptibility does not necessary 
coincide with the one in l|4.12|l and may, in particular, 
acquire an anomalous dimension 35] if d is low enough. 
However, this can only be due to infra-red divergent cor- 
rections to Eliashberg theory, which are fully captured 
by the effective low-energy model of Eq. (|4.13|) . The 
anomalous dimension emerges at the antifcrromagnetic 
QCP in d = 2 (Ref. but not in our case - 



Expression 


Definition 


Eq. 


Vf 


Fermi velocity 




m 


bare quasiparticle mass, m = I^f/vf 


12.211 


m 


effective (renormalized) quasiparticle mass 


EH 


TYlB 


band mass, determines the curvature of the Fermi surface 




g 


spin-fermion coupling constant 






ferromagnetic correlation length 




xoe 


value of the spin susceptibility at q — 




N 


number of fermionic flavors 




g=g 2 xo 


effective four-fermion interaction 


EH 


Nmg 

' TTV p 


Landau damping coefficient 


EH 


4:7TV p 


dimensionless coupling constant, it measures the mass enhance- 
ment: A = — — 1 

m 


EH 


373s 3 . . g 2 

8tt 3 iv~3, NE f 


frequency up to which E(w) dominates over iv in the fermionic 
propagator 


EH 




frequency up to which the fermionic and the bosonic mass-shells 
are well separated 


EH 


a- ~ g 

"fv'p NE p 


small parameter, measuring the slowness of the bosonic modes 
compared to the fermionic ones; the same small parameter jus- 
tifies the low-energy description 


EH 


r mN 


small parameter related to the curvature of the fermionic dis- 
persion 





TABLE I: List of the various parameters used in the text, their expression before the rescaling in N, and the reference equation 
where it is defined in the text. 



V. VALIDITY OF THE ELIASHBERG THEORY 

The essential part of the Eliashberg procedure is an 
a posteriori verification that the neglected terms in the 



self-energies are small. Quite generally, the validity of 
this procedure is based on the idea that the fermions are 
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fast excitations compared to the bosons, and hence the 
fermionic and bosonic mass shells are well separated in 
energy. 

When scattering off physical mass-shell bosons, the 
fcrmions are forced to vibrate on the bosonic mass shell, 
which is far away from their own. The electronic spec- 
tral function near the bosonic mass-shell is then small and 
this reduces the self-energy that arises from true fermion- 
boson scattering. In the case of the electron-phonon in- 
teraction, this is known as the Migdal theorem. 

The computation of the fermionic self-energy S(cj) 
gives us an idea of what the typical intermediate mo- 
menta and frequencies are in the problem. One can make 
sure that at criticality the typical fermionic momenta 



k — kp are of order Yi{w)/vf 



1/3 2/3 



jvp. On the 



other hand, the typical bosonic momenta q± along the 
direction of (i.e., transverse to the Fermi surface) are 
of the same order as the typical fermionic momenta, while 
the momenta qn transverse to kp (i.e., along the Fermi 



surface) are of order (7W) 1 / 3 » w t 



,V 3 „,2/3 







jvp. We see 

that for a given frequency u, the typical \q\ = \Jo\ + ?jj 

are much larger than k — kp, i.e. the effective boson ve- 
locity is much smaller than vp. One then expects that 
the Migdal theorem holds. 

The ratio of the typical fermionic k — kp and bosonic 
\q\ at the same frequency u is (ujQUj/^vp) 1 ^ 3 . At uj ~ Uq, 
this ratio then becomes: 



1/3 



9 



NEr, 



(5.1) 



and the slowness of the bosonic mode is then ensured in 
the quantum critical regime provided that ci« 1. This 
condition coincides, in our case, with the condition that 
the interaction should be smaller than the bandwidth. 
This is a necessary condition for the effective low-energy 
model to be valid, for if it is not satisfied, one cannot 
distinguish between contributions coming from low and 
from high energies. 

However, the smallness of a is not sufficient. In the 
integral for the fermionic self-energy, only one component 
of the bosonic momentum, q» , is much larger than k — kp, 
the other one is comparable: q± ~ k — kp. One needs to 
check whether the corrections to the Eliashberg theory 
scale as the ratio of k — kp to the modulus of \q\ or one 
of its components. 

To address these issues, we explicitly compute the ver- 
tex corrections and the fermionic self-energy at the two- 
loop level. 



Vertex corrections 



frequencies. We present the results below and discuss 
technical details of the calculations in Appendix C. 



1. Three-leg vertex with zero external momentum and 
frequency 

We begin with the simplest 3-leg vertex, with strictly 
zero incoming frequency fl and momentum q. The 
one-loop vertex renormalization diagram contains one 
bosonic line and is presented in Fig. [SJi. In analytic 
form, it writes: 



A3 
9 



q =n=o 
dujcPp 



'J dwd 2 p G(k F ,uj) 2 x(p,w) 



p 



P (if:(uj)-VFp t 



pi V 

2rn B J 



where we defined S(w) = to + T,(uj) and we have chosen 
kp along the x axis, so that p x = p± and p y = P\\- 

Since the poles coming from the fermionic Green's 
functions are in the same half plane, the integral over q x 
is only non-zero because of the branch cut in the bosonic 
propagator. At the branch cut p x ~ p v , so that we can 
safely drop the quadratic term in the fermionic propa- 
gators. Introducing polar coordinates, and integrating 
successively over the angle between kp and p, then over 
the modulus p, we obtain: 



9 



9 =O=0 



9 



ivp Jo 



du> 



S(w) 



(5.2) 



where w max is the frequency up to which bosons are slow 
modes compared to fermions, i.e. up to which £(w) jvp < 
(jlj) 1 / 3 . This frequency exceeds uo, so to find it we have 
to use E(w) w to. We then obtain 



y/NgEf 



(5.3) 



Note that for small values of a, w max 3> and the max- 
imum frequency up to which the bosons can be treated as 
slow modes well exceeds the upper limit of the quantum- 
critical behavior. 

Substituting S and w max into 1)5. 2[l . we obtain 
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(5.4) 



g=O=0 



This correction to the vertex can then be neglected pro- 
vided that a is small. 



We consider the vertex corrections due to the inser- 
tion of one bosonic propagator (three-leg vertex) and 
two bosonic propagators (four- leg vertex) . The behavior 
of these vertex corrections strongly depends on the in- 
terplay between the internal and external momenta and 



2. Three-leg vertex with finite external momentum 

We now turn to the three-leg vertex with zero external 
frequency but a finite bosonic momentum q. The one- 
loop renormalization is given in FigEl & n d its analytic 
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form writes: 



q,Q=0 



fl »/^pG(k F+P + q)W ) 

G(k F +p,w)x(p,w) 
duid 2 p 1 



l 



where p x is the projection of p along kp. 

As before, the integral over p x is can be reduced to 
the contribution from only the branch cut in the bosonic 
propagator. At the branch cut, p x ~ p y , hence we can 
neglect the quadratic term in the fermionic Green's func- 
tion, which then allows us to integrate over p y . Expand- 
ing in q x , and subtracting the constant term at zero mo- 
mentum calculated in (|5.2I) . we obtain (detail in the Ap- 
pendix C) 
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q,Q=0 



q=fl=0 



if/ ^log[zl]H] 

k F J\v Fqi \ M 

q x 



log \q x 



(5.5) 



When not only the bosonic momentum q is finite but 
also the external fermionic momentum for the static 
three-leg vertex is away from kp, the q x log dependence 
of the static Ag/g survives, but the argument of the log- 
arithm is the maximum of bosonic q and fermionic k — kp 
(we directed both along x). 

In general, the typical value of q x is much smaller than 
kp, hence the momentum dependent part of this vertex 
correction is small. However, we argue in the next sec- 
tions that because of the logarithmic term in (|5.5() . the 
insertion of this vertex correction into the SU (2) static 
susceptibility gives rise to a non-analytic q 3 ! 2 term. 

We also emphasize that although the calculations of 

look similar, the characteris- 



es 

9 



t;=S2=0 



and 

a 



q,n=o 



tic bosonic momenta are different for the two cases. At 
q = ft = 0, typical bosonic momenta in (|5.2|) are of or- 
der (7a;) 1 / 3 , i.e., near the bosonic mass shell. On the 
other hand, typical bosonic momenta in l|5.5|) are of or- 
der Ti(lo)/vf, i.e., near the fermionic mass shell. 



3. Generic three-leg vertex 

We next consider the same vertex with small but finite 
external momentum q and frequency Q. This diagram, 



presented in Fig. |2t>, reads 
A 5 



9 



/ dud 2 p G(k F + p + q, 

G(k F + p,w) x{P,v) 
dujd 2 p 1 



ft) 



V 



iE(aj) 
1 



VFPx 



iE(w + ft) - v F q x - v F p x 



Py 
2 m f 



,(5.6) 



where we have chosen ki? along the x axis, so that q x — 
qj_ and q y = q\\ . 

Integrating over p x first, one obtains two contributions: 
one arising from the poles in the fermionic Green's func- 
tions (which now can be in different half-planes since ft is 
finite), and the other from the branch cut in the bosonic 
propagator. The latter leads to the same result as (|5.2|l . 
up to small corrections due to the finiteness of the ex- 
ternal q and ft. Focusing on the other contribution, one 
has: 



A5 
9 



i— J <Lc j dp,. 



\Py\ 



+ \Py\ 3 



+ iE(w) - v F q x - aat 



(5.7) 



where the simplification of the frequency integral comes 
from the pole structure in p x . 

This generic correction to the vertex strongly depends 
on the interplay between the external q x , q y , and ft. 

When q = 0, but ft is finite, the integration over p y 
gives: 



Ag 
9 



U 



=o,n 



1/3 



did 



1 



VFl 



1/3, 



1/3 



E(fi 
= 0(1). 



■w) + £(w) 



(5.8) 



We see that the finite external frequency leads to a vertex 
correction that is not small, and one can make sure that 
higher-order corrections are of the same order. The series 
of vertex corrections have been summed up in |36| (for 
a charge vertex), where A.g/g was proved to diverge like 
(ujo/Q) 1 / 3 at low frequency, as expected from the Ward 
identity. 

When q and ft are both finite, 



A5 
9 



= T 



where 



VFq 

E(O) 



1 



, (1 - z) 2 / 3 + z 2 / 3 + ix 

has the following asymptotic behavior: 

f 1) = 0(1) 

1 ^»1) = OQ) 



(5.9) 



(5.10) 



(5.11) 
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If typical q is on the bosonic mass shell, then q 
and one has: 



( 7 fi) 



1/3 



Ag 
9 



q.Q 



E(Q) 



1/3 



(5.12) 



This is obviously small in a. 

It turns out that the behavior of the vertex correction 
is more complex and the result for Ag/g strongly de- 
pends on the direction of q compared to the direction of 
fermionic kp. This directional dependence is important 
for our purposes as we know from self-energy calculations 
that in the integral for the self-energy only y component 
of the bosonic momentum is near a bosonic mass shell 
and scales as (7SI) 1 / 3 , the x component of bosonic mo- 
mentum is actually of the order of H(u>)/vf, he., is near 
a fermionic mass shell and is much smaller. We therefore 
take a more careful look at the vertex correction. 

For the case where q x ~ E(f2)/«F and q y ~ ^fi) 1 / 3 , 
one would argue from (|5.12|) that the vertex correction 
now becomes of order 0(1) and is no longer parametri- 
cally small. However, the computation that lead to (|5.12|) 
cannot be extended to the strongly anisotropic case as for 
external q x ~ E(f2)/«p and q y ~ (-/fl) 1 ^ 3 , the curvature 
of the fermionic dispersion becomes relevant and changes 
the result. The full dependence on q x and q y is rather 
complex and we restrict ourselves to the case when 







1 m B 
N m 



< 1. 



(5.13) 



the 



In this situation, q 2 /rriB ~ (vpq x )/l3 3> VFq x , he 
the quadratic term in fermionic propagator dominates 
Performing the integration, we then find that: 
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fl 2 



qS2 



7f7 



2/3 



log 2 



( 7 0) 



1/3 



/? 2 log 2 /3, 



(5.14) 



It follows that even when only one component of the 
bosonic momentum is near a bosonic mass shell, the ver- 
tex correction is small if (3 is small. This is the second 
condition for the Eliashberg theory to be controllable at 
criticality. 

The smallness of j3 can be ensured by either extending 
the theory to large values of N, or by considering a very 
strong curvature of the Fermi surface which implies that 
tub *C m. Even though the latter can hardly be satisfied 
for realistic Fermi surfaces, we emphasize that the curva- 
ture of the dispersion plays a crucial role in the theory, 
for even in the case of N 3> 1, the vertex correction is of 
order O(l) without this curvature. 




k F ,0 



k F ,0 



k F ,0 



FIG. 2: Three-leg vertices: a) zero external momentum and 
frequency b) finite momentum c) generic vertex. 



the curvature of the Fermi surface. This leads to a vertex 
of order O(l) even in the large-N limit and the pairing 
problem then has to be carried out exactly within the 
Eliashberg theory. 

This vertex renormalization is presented in Fig. [2 and 
its analytic form is given by: 



A5 
9 



Cooper 



1 J dud 2 q x(q, to) G(k F + q, to) 

ft) 



xG(-k F 
dojePq 



q, — to - 
1 



1 



q 2 i£(w) - v F q x - 
1 



-iE(w + Q)- v F q x - 



(5.15) 



Integrating over q x , restricting ourselves to the con- 
tribution from the fermionic poles (the one from the 
branch cut can be proved to be smaller), we find that 
the quadratic terms cancel out, leaving us with: 



Ag 

g 



D 



dq y q y 



Cooper V F J\a\ E(W + 0) + E(W) Jo 1^ + 1y ' 

(5.16) 

Performing the remaining integral, the pre-factor sim- 
plifies and we obtain: 



A5 

g 



9 



Cooper 



1/3 V 3 

n 

D 



log 



log 



(5.17) 



where we assumed that we were in the quantum critical 
regime, i.e. |0| < c^o- 

We emphasize that the pre-factor of the log in (|5.17|l 
is O(l), even when one takes into account the curvature 
of the fermionic dispersion. The result of Eq. (|5.17() 
confirms previous studies :Sj advocating that the system 
at a ferromagnetic QCP is unstable towards pairing. 



Four-leg vertex 



4- Pairing vertex 

By contrast to the previous vertices we studied, the 
pairing vertex in the Cooper channel is not sensitive to 



We consider now higher-order corrections to the vertex 
through the example of a four-leg vertex correction with 
two crossed bosonic lines, also called a Cooperon inser- 
tion. Analytically, the expression for this renormalized 
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a) 



o,-a 



k F ,0 



>k F + q, Q 



b) 



B. Self-energy corrections 

1. Corrections to the self-energy at the two-loop level 



kp,0 k F + q,Q k F ,0 

FIG. 3: a) Cooper pairing vertex b) Four-leg vertex 

vertex, presented diagrammatically in Fig. |3J writes: 

'Q + uj p + q N 



r 2 (g,fi) ~ g 2 J dujd 2 p x s 



2 1 2 



'0,-uj q-p\ fCl + w p + q 



xG['^k F ■ q P 



2 ' 2 

' (5.18) 
After integrating over p x (projection of p along k F ), and 
uj, were arc left with: 



r 2 (g,fi) 



.9 2 n 



dz 



(- + «)' 



- 4z 2 



«f|«vl 3 ./ 2iE (§) - i^fe - &(1 



1 



(z 2 + 2z - 



3/2 



( z2 ~ 2z + %) 



3/2 



7 |0| 

(5.19) 

where we have changed p y into z = Py/\<ly\- 

This renormalized 4- leg vertex T2(q, f2) has to be com- 
pared with the bare four-leg vertex, whose analytic form 
is given by the bosonic propagator multiplied by g 2 : 



ri(«,n) 



fj 



2l£i' 

9 



(5.20) 



In the case of a typical external momentum q x ^ q y ^ 
(7O) 1 / 3 , the ratio T2/T1 is of order a. For a typical ~ 
Ti(oj)/vf and q a ~ (7W) 1 / 3 , we obtain, to logarithmic 
accuracy, 



We found in our analysis of the vertex corrections 
that the result depends on the interplay between the 
typical momentum and frequency. In our estimates, 
we considered two regions of external q and tt, namely 
Qx ~ q y ~ (7^) 1/3 and q x ~ t(uj)/v F , q y ~ (7ft) 1 / 3 . In 
both cases, we found that vertex corrections are small. 
We verify here that the two-loop self-energy, obtained 
by inserting vertex corrections into one-loop self-energy 
diagram, is also small. 

The two-loop self-energy diagram is presented in fig. 
□ We have: 



£ 2 M 



J dujid 2 qi y dw 2 d 2 g2X(qi, ^i)x(q2,^2) 

x G(k F + qi, uj + wi)G(k F + q2, ui + uj 2 ) 
x G(k F + q! + q 2 , uj + uo 1 + uj 2 ) 

qi 92 



g 2 I dujid 2 q\ j duj2d 2 q2 



1 



7|wiJ +qf j\uj 2 \ + q 2 
1 



*E(w + UJ!) - v F q lx + W2 ) _ Vpq2x _ ^_ 

1 



iT,(uj + UJ1+ uj 2 ) ~ v F q lx - v F q 2x - ^ 2m 



>.m B 

(5.22) 



where we recall £(w) = uj + E(w). 

Integrating successively over and g 2x , and rescal- 
ing the remaining momentum components by introducing 
x = <?iy/ (tI^i l) 1/3 and V = 92y/(7|w2|) 1/3 , we obtain: 



£ 2 (w) 



m B g 



duj 2 



U F Jo 
dxdy 



duj\ 



(7 2 wiw 2 ) 2 / 3 



(5.23) 



xy + zC(l + kl 3 )(l + |y| 3 )' 



( 7 2 ^ 2 ) 1/3 



where C = 777,5 

As the typical frequencies ui\ and o> 2 are of order uj, 
the typical value of Q is of order /3 -C 1. Expanding then 
in l|5.23|l to first order in £ and performing the remaining 
integrals, we obtain in the quantum-critical regime: 



r 2 

J- 1 



(5.21) 



This last result again critically depends on the curvature 
of the Fermi surface: neglecting the quadratic terms in 
the fermionic propagators, one would obtain r 2 /ri = 
O(l). We see that likewise to the three-loop vertices, the 
smallness of the crossed vertex T 2 (g, fl) requires both a 
and to be small. 



E 2 (w) 



mg 2 r r ciog 2 c 

32 . „2 



E(w)^logr«/3, 

,V 3 „ 2 / 3 



(5.24) 



o; ' W' " is the self-energy in the 



where E(w) = Ei(tj) 
Eliashberg theory. 

This result agrees with the one obtained in 2 



shows that E 2 (w) ~ Ei(w) x 



As 



q,Sl 



where — 



and 
is 



q.n 
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given by i|5.14[) . This last result implies that typical inter- 
nal q and f2 for the Eliashberg self-energy and for T, 2 (uj) 
are the same. 

It is also instructive to compare these two-loop re- 
sults, obtained as an expansion around the Eliashberg 
solution, to the perturbation expansion around free 
fermions. In the latter, we found in Eq. I|3.12[l that 
E2(w) ~ uj log 2 10 whereas in the former we have T,2(uj) oc 
j3(ul /3 oj 2 / 3 ) log 2 (3, Eq. The free-fermion re- 

sult can be reproduced if we neglect the self-energy in 
(|5.23l) . We see that the expansion around free fermions 
does not reproduce the correct frequency dependence of 
£2(0;). This obviously implies that if one expands around 
free fermions, there exist higher-order terms associated 
with insertions of the self-energy into the internal 

fermionic lines, which may overshadow the two- loop re- 
sult around free fermions. Accordingly, near the QCP, 
the expansion around free fermions does not converge, 
even if the curvature of the fermionic dispersion is in- 
cluded. On the other hand, the expansion around the 
Eliashberg solution is regular and holds in powers of the 
small parameters a and (3. 



2. Momentum dependence of the self-energy and the 
density of states 

Along with the vertex corrections, we also neglected 
the momentum dependence of the fermionic self-energy 
in order to proceed with the Eliashberg scheme. We now 
verify whether the momentum dependent part of the self- 
energy £(k, u) = 0) = E(k) remains small when evaluated 
with the full fermionic propagator. The k dependent self- 
energy is given by: 



E(k,0) = 3ig 



Kg 2 J ^^G(k + q,fi) X (q,0) 
3ig f dflcPq q 



(2tt) 3 J iE(fi) - e k+q 7 |n| + q 3 



. (5.25) 



Defining the angular variable as e k +q = e/c + vpqcosO, 
integrating over it, and expanding to linear order in ek 
we obtain: 



£(k,0) = -3ie k x 



/ 



2?r- jo 
q 2 dq 



dfl E(fi) X 



( 9 3 + ( 7 |fi|))(W?) 2 + £(n) : 



3/2 ' 



(5.26) 



A simple experimentation with the integrals shows 
that the integration over momentum is confined to q ~ 
£(f2)/wir, while the frequency integral is confined to 
Q < ujmax, where Lu max , defined in 1|5.3|1 . is the scale 
where (jQ) 1 / 3 = E(fi)/uj?. The computation of E(k) is 
given in the Appendix B, and the result is 



E(k,0) 



-"/■ x 3 ^ X i' 3 o 3( ! 8 V" = -itk 0.253Va (5.27) 



Using our definition of E, we then obtain that E(k, 0) 
gives rise to a small, regular correction to the quasipar- 
ticle mass e k — i£(k, 0) = e| = v F (k — kp). 



v F 



(1 - 0.253y/a) 



(5.28) 



Like the vertex correction at zero external bosonic mo- 
mentum and frequency, this small correction is of order 
0(y/a) and comes from frequencies of order a; max . 

The momentum dependent self-energy, unlike the fre- 
quency dependent part, generally gives rise to corrections 
to the fermionic density of states (DOS) 

de k 



N(w) 



-N 



-ImG(e k ,u)) 



(5.29) 



where iV is the DOS of free fermions. 

Assuming that E(fc) is small and expanding in it in the 
fermionic propagator, we obtain, in real frequencies 



N(u) 
No 



= 1 — Im 



E(k) 



(5.30) 



t£(- 



Sunstituting (|5.27|) into l|5.3UI we find that the density 
of states just shifts by a constant. In order to extract 
the frequency dependence of the density of states, one 
has to evaluate the momentum-dependent self-energy to 
next order in e k and on the mass shell, where e k = iE(w). 

The evaluation of the self-energy near a mass shell 
generally requires extra caution as the self-energy may 
possess mass-shell singularities j^S). We, however, have 
checked in Appendix [0 that in our case the self-energy 
does not possess any mass-shell singularity, and the self- 
energy remains finite on the mass shell. 

The calculation of the self-energy to order t\ is dis- 
played in the Appendix B, and the result is 

E(k,.)^x^Mlogfi (5.31) 
87rA hjp 

Substituting this self-energy into the expression for the 
DOS and converting to real frequencies, we find at small 

N(u) ~ N (l - (£) ^ log , (5.32) 



where we explicitly defined 



128tt 5 / 2 



N 2 E 2 F 



(5.33) 



2V2tt 3 / 2 



0.45 3 / 2 3 3 / 4 g ' 

The correction to the fermionic DOS was earlier com- 
puted by Lawler et al 01 using the bosonization tech- 
nique. They obtained the same uj 2 ^ 3 logo; dependence 
as in H5.32I The agreement with Ref. [lj] is, however, 
likely accidental as they evaluated N(lj) in the expan- 
sion around free fermions, while we obtained N(u>) by 
evaluating the k— dependent self-energy at e k = iH(-iui). 
If we neglected the Eliashberg self-energy (i.e., expanded 
around free fermions), we would obtain u> logw correction 
to the DOS. This last result agrees with the one obtained 
by plj using the same technique as in [T^ . 
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C. Summary 

We have shown in this section that there are two con- 
ditions for the validity of the Eliashberg theory that one 
can recast as the smallness of two parameters: 



(5.34) 

The first condition is quite generic for a low-energy 
theory since it requires that the fermion-fermion inter- 
action mediated by the exchange of a boson should be 
smaller than the Fermi energy. Otherwise, the physics 
is not restricted to the vicinity of the Fermi surface any- 
more. The parameter a plays the same role as the Migdal 
parameter for the electron-phonon interaction: it sets the 
condition that fermions are fast excitations compared to 
bosons. In the scattering processes that are small in a, 
fermions are forced by the interaction to vibrate at fre- 
quencies near the bosonic mass shell. They are then far 
from their own resonance and thus have a small spectral 
weight. 

However, the condition a <C 1 is not sufficient to 
construct a fully controllable perturbation expansion 
around the non-Fermi liquid state at the QCP. In spa- 
tially isotropic systems there exist vertex corrections for 
which the external momentum has a component on the 
fermionic mass-shell. These corrections don't contain a. 
However, these corrections are sensitive to the curvature 
of the Fermi surface, and are small if (3 is small which can 
be achieved either by imposing tub < m or by extending 
the theory to a large number N of fermionic flavors. 

A word of caution. In evaluating the rcnormalization 
of the static vertex, we silently assumed that y/a <C /3, 
i.e., g/Ep < (mB/ni) 2 /N. At very large N, this is no 
longer valid. For this situation, i.e., when (3 <C \fot, our 
estimates show that the static vertex is even smaller than 

\[0L. 

Finally, the pairing vertex in the Cooper channel stays 
of order O(l), signaling the possibility of a pairing insta- 
bility close to the quantum critical point. Nevertheless, 
we assume, based on explicit calculations worked out in 
, that the quantum critical behavior extends in the pa- 
rameter space to a region where the superconductivity is 
not present. 



VI. INSTABILITY OF THE FERROMAGNETIC 
QUANTUM CRITICAL POINT 

We found that the Eliashberg theory for fermions in- 
teracting with gapless long-wavelength bosons is stable 
and controlled by two small parameters. We verified this 
by calculating the fermionic self-energy in a two-loop ex- 
pansion around the Eliashberg solution 

One may wonder whether the same conclusions hold 
for the bosonic self-energy as well. In particular, what 
are the corrections to the static susceptibility Xs(q, 0)? 



Naively one could assume that they are unimportant and 
do not change the bare q 2 behavior of the inverse bosonic 
propagator at the QCP. 

For a ferromagnetic SU(2) QCP, for which the mass- 
less bosons are spin fluctuations, we show in this section 
that the corrections to the static spin susceptibility are 
non-analytic: they scale like q 3 ^ 2 , and do not contain any 
pre-factor except for a proper power of kp. Such terms 
obviously overshadow the regular q 2 of the bare suscep- 
tibility at small enough momenta. These terms therefore 
belong to Eliashberg theory, which has to be extended to 
incorporate them. 

The physics behind the q 3 ^ 2 term in x(<Zi 0) at a 
ferromagnetic QCP is, by itself, not directly related 
to criticality: far away from the QCP, the spin sus- 
ceptibility also contains negative, non-analytic \q\ term 
El El El HH This term gradually evolves as 
the correlation length £ increases, and transforms into 
the q 3 / 2 term at the QCP. Both these non-analyticities, 
at and away from the QCP, emerge because the boson- 
mediated interaction between fermions contain a long- 
range dynamic component, generated from the Landau 
damping. 

For charge fluctuations, the q 3 ! 2 terms appear in the 
individual diagrams for the susceptibility but cancel out 
in the full formula for x(<?, We discuss the physical 
origin of the difference between spin and charge suscep- 
tibilities in the next section. 

One of the reasons why the non-analyticity of the static 
spin susceptibility at the QCP has not been analyzed 
much earlier is because it was widely believed that an 
itinerant fermionic system near a ferromagnetic QCP is 
adequately described by a phenomcnological 2 + ID 4> 4 
bosonic theory (in our case, the role of 4> ls played by 
the vector field S) with the dynamic exponent z = 3 and 
a constant pre-factor for the (f> term lj. In dimensions 
d > 4 — z = 1, the model lies above its upper critical 
dimension and the (j) 4 term is simply irrelevant. 

In this section, we derive the effective </> 4 theory from 
the spin-fcrmion model Hamiltonian, and show that it 
contains two new elements absent from the phenomeno- 
logical approach. First, the pre-factor of the (j) 4 term 
strongly depends on the ratio between the external mo- 
menta and frequencies, and contains a non-analytic term 
in addition to the constant one. Second, there also exists 
a cubic 4> 3 term whose pre-factor, although vanishing in 
the static limit, also strongly depends on the interplay 
between the external momenta and frequencies. We can 
recast the non-analytic q 3 ! 2 term in the static spin sus- 
ceptibility as arising from these cubic and quartic terms 
in <p. 

We also prove that the non-analyticity appears in 
the temperature-dependent uniform static susceptibility 
Xs(T). e show below that xj 1 ^) oc T\ logT|, again with 
a negative pre-factor. 

Finally, we show that the instability of the ferromag- 
netic QCP can also be seen from the fermionic self-energy, 
which acquires singular terms beginning at the three-loop 
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0,0 



FIG. 4: "</i> 4 " and V type diagram 



order in the case of spin fluctuations. 



A. Hertz-Millis-Moriya theory revisited 

In order to derive a quantum critical cj) 4 model, one has 
to integrate the fermions out of the partition function, 
noticing that the Lagrangian of the spin-fermion model 
is quadratic in the fermions. 

Expanding then in the number of bosonic fields S, the 
quartic term in the effective action reads: 



/ 



d 2 qd 2 pd 2 p' dVLdvdv 1 



^(p,p',q, v, v',Q) 



+ &p+q/2 ' ^>-p'-q/2 Sq/2-p ' ^p'-q/2 



(6.1) 



where p,p' ,q are bosonic momenta, and v, z/,f2 are 
bosonic frequencies. 

Our goal here is to prove that the pre-factor A is not a 
regular function of momenta and frequencies. To simplify 
the presentation, we choose to study only the dependence 
of A on q and fl and set p, p', v, v' to zero (see Fig0J). 

The analytic form of this prefactor then writes: 



A(q,Q) ~ Ng* J dcu J d 2 k G(k,w) 5 
x G(k + q,w + fi) 2 . 



(6.2) 



Defining 9 as the angle between k and q, and perform- 
ing the integration over £fc,and the angular variable, we 
find: 



A(q,il) 



Nmg 4 

LUqQ 



dz 



(^) 2 -2(z 2 /3 + (1 _ z)2 / 3 y 



(^) 2 + (- 2/3 + ( 1 --) 2/3 )' 



5/2 



, (6.3) 



where we defined z — ui/fl, and neglected at this stage a 
regular part that comes from large values of efe and for 
which the curvature is relevant. 

We see that A(q, SI) depends on the interplay between 
the momentum and frequency. We can identify two 
regimes 

• If \q\ ~ (7IJ1I) 1 / 3 , i.e. if the bosonic momenta are 
near the bosonic mass shell, the self-energy in the 
denominator can be neglected, the frequency fac- 
tors in the numerator and the denominator are can- 
celed out, and we obtain: 



A(n) 



Nmg 4 



1 



(6.4) 



We see that A can be safely replaced by a con- 
stant OL a )- This is consistent with the previous 
works pj. The agreement is not surprising as the 
relation q ~ (7|f2|) 1 / 3 is assumed in power count- 
ing based on z — 3 scaling. Note that the condi- 
tion \q\ ~ (7IS1I) 1 / 3 only specifies the magnitude 
of q, one of its components (e.g., q x ) can be much 
smaller. 

If q ~ S(w) ~ lu^ 3 Q 2 / 3 vf (at oj < ujq), i.e. when 
a boson resonates near a fermionic mass shell, the 
z = 3 scaling arguments are not applicable. We 
have in this regime 



A(q,Sl) 



Nmg 4 Nm 1 w n 



w f2 xl Va ft 



(6.5) 



In this case, A(fi) is a singular function of fre- 
quency, and cannot be replaced by a constant. 

We see therefore that the pre-factor of the 4> 4 term is 
actually singular outside the scaling regime of a z = 3 
theory. 

In the similar spirit, one can construct a cubic term in 
the bosonic fields. 



d 2 qd 2 pdQdv 
J2nf 



B(p,q,u,fl) S p ■ (S g _ p/2 x S_ g _ p/2 ), 

(6-6) 

where the pre-factor B is a convolution of three fermionic 
Green's functions as presented in Fig^ and is given by: 

B(q,p,n,i/) ~ Ng 3 J dw J d 2 k G(k-p,u-v) 
xG(k-p/2-q,w-i//2-n) G(k,w). (6.7) 

Proceeding as for the quartic term, we set for simplicity 
p = 0, v = 0, and integrate over e& and the angular 
variable, leading to: 



Nmg 3 



2/3 



dz 



{z 2 ' 3 + (1 - z) 2 ' 3 ) 



3/2 



(6.8) 
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where we again introduced the rescaled frequency z — 
ui/Q. 

We can again identify two regimes. 

• In the z = 3 regime where q ~ ^fi) 1 / 3 , one can 
expand for large -jj^fy, which leads to: 



u( n\ Nm 9 ( n 
B{q, 0) a 



Xo 



2/3 



(6.9) 



where io max is given by (|5.3|) . This term is small 
in the quantum critical regime where f2 < wo ~ 
w, MI a 3 / 2 , and can be safely neglected. 

• For q ~ S(fi), the remaining integral is of order 
0(1) and the result writes: 



which is large and cannot be neglected. 



(6.10) 



We demonstrate in the next section how the singular 
behavior of the cf> 3 and 4> 4 terms leads to a non-analytic 
contribution to the static spin susceptibility. 



B. Non-analyticity in the static spin susceptibility 



We now estimate the effect of these singularities on 
physical quantities. Both the (f> 3 and <p terms in the 
effective action give rise to corrections to the <fi 2 term, i.e. 
to the spin susceptibility. These corrections are obtained 
diagrammatically by contracting the external legs of the 
4> 3 and </> 4 terms, as shown in Fig. The computations 
are described in detail Appendix IdI 



The contributions from the (f> 4 terms have been con- 
sidered in our short publication |3lj| . The contributions 
from cubic terms were missed, and were first considered 
m in the analysis of the spin susceptibility in the 

paramagnetic phase, away from a QCP. 



In analytic form, we have, using x(q, 0) = Xo/(£ 



U(q, 0)), and II = U^q, 0) + U 2 {q, 0) + n 3 (g, 0) 



n 4 (g,0): 



III (g,0) 

n 2 (ff,o) 
n 3 (g,o) 



it 



= r 



= r 



Nf 
(2^) 6 Xo 
s Nf 

(2tt)6xo 
s N 2 g 3 

3 Wfxl 

G(oj + Q,k + l + q) G(uf, k') G(wf, k' + q) G{u' + n,k' + l + q) 



d 2 kdujd 2 ldfl X s(l, fi) G(u, k) G(u + Q,k + l) G(u + il,k + q + l) G(u, k + q) 



d 2 kdcud 2 ldn Xs (l, G(oj, k) 2 G(uj + O, k + I) G(u>, k + q) 



d 2 kdcud 2 k'duj'dldn Xs(l, O) Xs(q + I, O) G(w, k) G(lu, k + q) 



n 4 ( g ,o) = r 



N 



2^3 



d 2 kdujd 2 k'duj'dldn Xs (l, Xs(l, 0) G(oj, k) G(ui, k + q) 



4 i^fxl 

G(uj + n,k + l + q) G(oj', k') G(oj' + n, k' + I) G(w' + Sl,k' + l + q) 

I 



(6.11) 



The factors of N come from the fermionic loops and 
the numerical pre-factors from the following spin sum- 



mations: 



La — 



a a pap 7 a lS as a = -2 
^2 CT "/3 (°>y ■ ^Se) {<r-ta ■ cr c5 ) crf c = 8 

a,/3,7,<5,e,£ 



rf = 2J a af3 (°>y " Vet) (o" 7 a ' ^5c) a fe = _8 

(6 - 12) 

These four diagrams are related by pairs. To verify 
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a) 9 '^JJh b ) 9 '° 

c) ^'^{^j ^ ^''^^^^^^ 

FIG. 5: Corrections to the polarization bubble from diagrams 
with one and two extra bosonic lines 




this, it is useful to expand the products of Green's func- 
tions according to: 



G(a;i,fei)G(o; 2 ,fe 2 ) 



. (6 : 13) 

Applying this to II2 (g, 0) , we find that it splits into 
two parts. In one part, the poles in tk are located in 
the same half-plane, leading to a vanishing contribution. 
The remaining term in LI2 (g, 0) is related to IIi(<7, 0) in 
such a way that: 

ll 1 (q,0) = — =±ll 2 (q,0). (6.14) 
I 2 

(see Appendix IdI for details). 

Similarly, 113(9, 0) an d 114(9, 0) are related as 

n 3 (g,o) =-^114(9,0). (6.15) 
1 4 

Collecting all four contributions and using the relations 
between pre-factors, we obtain: 



U{q,0) = Tl A (q,0)+Tl B (q,0) , 

n A ( q ,o) = n 1 ( <z ,o) + 2n 2 (< z ,o) 

Nf 



(6.16) 
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(2tt) 6 



d 2 Kdujd 2 ldtt Xs(l,tyG{uj,k)' 2 



xG(oj + n,k + 1)G(uj, k + q) 
U B (q,0) = n 3 (q,0)+n 4 ((7,0) 
N 2 f 

= 16 



. , d 2 kdud z k'dLj'didn y s a,o) 

x X s(q + l, 0) G(lu, k) G(uj, k + q) 
xG(Lu + n,k + l + q) GV, k' + q) 
xG(w', k') G(oj' + n,k' + l + q) (6.17) 



1. Fermi-liquid regime 

Away from criticality, the correlation length £ is finite, 
and at low frequency, the system is in the Fermi-liquid 
regime. The fermionic self-energy is S(w) = Xuj, Eq. 
(IP) . 



The spin susceptibility in this regime has been ana- 
lyzed in [H IHlH IHlH H3 It was shown there that to 
the lowest order in the interaction, ils(g, 0) = Ha(<7, 0), 
i.e., Il(<?, 0) = 211,4 (?,0). Beyond leading order, ns(g,0) 
and H A (q, 0) are not equivalent but are of the same sign 
and of comparable magnitude. For simplicity, we assume 
that the relation Hsiq, 0) = 11,4 (q, 0) holds in the whole 
Fermi liquid regime. We will see below that even at the 
QCP, 11^(5,0) and U.A(q, 0) are quite similar (at QCP 

n B (g,o)«i.3n A (g,o)). 

Introducing cos 9 = -^p and cos6>' = rgA , and suc- 
cessively integrating over |k|, ui and 9', i|6.16[l can be 
rewritten as: 



n(g,o) 



Q 



TT 3 (l + X)v F J 



dz 



dO 



cos sin i 



1 



tan 4> {i sin <j> ■ 



cos a cos ( 



(6.18) 



y/l + z 2 (s'm(j> + i cos <fi cos 9) 2 ' 

where we defined 7 = y^j , and introduced the new vari- 
ables z and cj), which satisfy z cos <fi — ^ and z sin <p = 

(1+A)Q 

The universal part of LT(g, 0) can be isolated by sub- 
tracting from it the constant part 11(0,0). The integral 
over z then becomes convergent. Integrating successively 
over z, 6 and 9, we obtain: 



n(<?,o) = - 



fj 



ir 2 v F (l + A) 



\q\ H 



1 + A 



(6.19) 



where H(0) = §, and H{x > 1) » 2/(3x 2 ) We do re- 
cover the non- analytic \q\ correction to the static spin 
susceptibility in D = 2, as obtained in earlier studies 

Note that Eq. 
liquid relation Xsil 



does not contradict the Fermi 
0,u = 0) cx (1 + F M )/(1 + F , Q ), 
where iq. s and Fo. a are Landau parameters. The Fermi 
liquid theory only implies that the static spin suscepti- 
bility saturates to a constant value as q — ► 0, but does 
not impose any formal constraint on the q— dependence 

ofx s (g,o>). 

As one gets closer to the QCP, A = 3g/(47ru5^ 1 ) di- 
verges and the pre-factor of the \q\ term vanishes as: 



n(g,o) 



16 
9^ 



r x i?i 



(6.20) 



This is not surprising since the Fermi liquid regime ex- 
tends on a region of the phase diagram that shrinks and 
ultimately vanishes as one approaches the QCP. 
Now, two different scenarios are possible: 

• the divergence of £ at the QCP completely elim- 
inates the non-analyticity and the expansion of 
II(<7, 0) begins as q 2 , like in a bare spin suscepti- 
bility, 
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• the self-energy T,(ui) oc w 2 / 3 at the QCP still leads 
to a non-analytic term II(g,0) oc \q\ s , with 1 < 6 < 
2, which dominates over the bare q 2 term. 

We show in the next subsection that the second scenario 
is realized, and at the QCP, one has IT(q, 0) oc Igl 3 / 2 . 



2. At criticality 

At the QCP, we have to take into account two new 
elements: the bosonic propagator is massless (£,~ 1 = 0) 
and the fermionic self-energy is no longer Fermi-liquid- 
like, it is given by (|4.9|) . 

The full calculation of II^^O) and Hsiq, 0) is pre- 
sented in Appendix [Dj We just outline here the main 
steps of the computation and show where the q 3 ^ 2 term 
comes from. 

Consider first Ha(<1i 0). Using l|6.16[) as a starting 
point, and substituting the full form of the spin suscep- 
tibility, Eq. i|A8() . we then expand £k+i and £k+ q , and 
integrate successively over l y (projection of 1 perpendic- 
ular to kp) and e^, leading to: 



U A (q,0) = lGi 



mg 



271 



d6 



+oo 



dL 



-\-oo 



d.9. 



du 



-Q 



^g + c 2 E 2 (») 

( 7 |n|)V3 



X (7|^|) 1 / 3 i£(w +0) - - v F l x + v F qcos6' 

(6.21) 

where c ~ 1.19878 (see IA6|) . and h(x) is the spin suscep- 
tibility X(Z,fi)/xo integrated over the momentum com- 
ponent l y . The asymptotic behavior of h(x) is given by: 



h(x » 1) = I 

h{x < 1) = ^ + (log 2 - l)x 2 - 



x log X 
2 



(6.22) 



Since the integrand in l|6.21|l has poles in l x located in 
the same half-plane, the only non- vanishing contributions 
to II comes from the non-analyticities in h(x). 

There are two non-analyticities in h(x). The first one 
comes from the l/x behavior at large x, which extends to 
x ~ 1. This is a conventional non-analyticity associated 
with bosons vibrating near their own mass shell, since at 

X ~ 1, l X ~ ly~ (7|"|) 1/3 . 

Subtracting the universal U^a (0, 0), expanding in q in 
(|fT2T|) and substituting l x ~ (7IOI) 1 / 3 , we find for this 
contribution to II: 



(a) 



oc q 



nig 

u F 7 5 / 3 



dQ 
ft 2 / 3 



oc 



Vaq 2 



(6.23) 




FIG. 6: Integration contour: the hatched region stands for a 
branch cut, and the cross for a pole. 



where ft ma 

We see that the integration over the momentum range 
relevant to the z — 3 scaling regime yields a regular q 2 
contribution to the static susceptibility. Not only this 
contribution is regular in q, but it is also small in a. 
This result is similar to the one we obtained in H5.2(l for 
the static vertex at a vanishing momentum. 

However, Eq. H6.22J) shows that h(x) has a non- 
analytic x 2 log 2; term already at small x, i.e. far from the 
bosonic mass shell: the branch cut associated with the 
logarithmic term emerges at vfI x ~ The typical 

value of l y in the integral that leads to this x 2 log x term 
is also of the same order, although larger in the logarith- 
mic sense. This implies that this second non-analyticity 
comes from a process in which the bosons are vibrating 
near a fermionic mass shell and far from their own. 



Furthermore, this logarithmic term in 16.22f) comes ex- 
clusively from the Landau damping term in the bosonic 
propagator - the q 2 term in x(h ^) can be safely omitted. 
Indeed, one has: 



dL 



-A 



7l"l 



/(A) 



2 7 ft 



(6.24) 



Substituting the logarithmic term from (|6.22l) into 
(|6.21(l . we subtract the non- universal constant term 

11^ (0,0), which makes the integral over l x convergent. 
Introducing z = / x /(cS(0)), one can perform the inte- 
gral over z over the contour of Fig. which leads to the 
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following contribution to II: 



2 r+oo 



ni%,o)-<>(o,o) = i^ 2 / del dy 



dw- 



cos 2 8 



S(fi) 2 7 [ c -l ((! _ u,J2/3 +w 2/3) +y ] 


[y + c-l ((1 - «0*/3 + ^/3)] 2 + (^ii) 

(6.25) 



where we defined w — w/fl. 



3/2 



and 



Introducing; the new variables t = ( 77 

\VFqcos9 

v = t [y + cT l ((1 — w) 2 / 3 + u> 2 / 3 )] , it becomes possible 
to separate the integrals, leading to the following final 
result: 



i4%,o)-n^(o,o) 



(6), 



-0.8377? 



3/2_ 



Nmg 2 



4 3/2 1/2 



= -0.1053 V fc F q 



3/2 



(6.26) 



We emphasize that this dependence comes from bosonic 
modes vibrating at the fermionic mass-shell. This ex- 
plains why the result of (|6.26|l is not small in a, as this 
small parameter measures the softness of the mass-shell 
bosons compared to the mass-shell fermions. 

The integrals for Hb(q, 0) cannot be exactly evaluated 
analytically, but an approximate calculation in Appendix 
E yields 



l4%,0) -11^(0,0) = -O.Mg 



r(6) 



3/2 



PF, 



such that the total 



U(q, 0) - n(0, 0) = -0.25g 3/2 v^F 



(6.27) 



(6.28) 



We see that H(q, 0) is still non-analytic in q and the 
prefactor is negative. At small q, the negative q 3 ! 2 term 
well exceeds the regular q 2 term. This implies that the 
static spin susceptibility is negative at small momenta, 
i.e. a ferromagnetic QCP is unstable. We discuss the 
consequences of this instability in the concluding section. 
The momentum q m i n below which \s is negative is de- 
termined by 



X 



cx q 



^■ 2 ^LVpf 



0, 



(6.29) 



which gives q m in — 0.06fcp. Parametrically, q m in is 
of order kp, which is the largest momentum scale in 
our problem. Strictly speaking, this suggests that the 
whole quantum-critical theory for the ferromagnetic case 
is not valid, for the quantum critical behavior extends 
up to energies of order ojq, i.e. up to momenta of order 
q < ^o/vf ~ a 2 hp <C kp. Numerically, however, q m in is 
much smaller than kp. This implies that for reasonable 
values of a, there exists an intermediate momentum and 
energy range q m in < q < ^o/vf where the system is in 



the quantum-critical regime, but the static spin suscep- 
tibility is still positive. 

The q 3 / 2 non-analyticity can also be viewed as emerg- 
ing from q log q momentum dependence of the static ver- 
tex. Using Eq. (|5.5|l . one can rewrite: 



II(g,0) ~ Ng / d 2 kdio 



9 



iS(w) - e fc iS(w) - e k+q 



(6.30) 



Performing the contour integration over e^, and changing 
variables into y = efe/Eu; and t — ^Jluquj j '((vpq) 3 / 2 ), we 
obtain: 



U(q,0) ~ yjkpq 



,3/2 



(6.31) 



We note in this regard that the non-analytic momentum 
dependence of the fermion-boson static vertex also comes 
from bosons vibrating near fermionic mass shell, i.e., it 
emerges due to the same physics as we outlined above. 



C. Temperature dependence of the static uniform 
spin susceptibility 

In this section, we show that the static uniform sus- 
ceptibility is negative at finite temperature above a fer- 
romagnetic QCP. To demonstrate this, we compute the 
static uniform Ua(<7 = 0, u> = 0, T) = n^(T), assuming 
that = 0. The contribution from 11b (T) is of the 
same sign and comparable in magnitude. We have 



n.iiD.D Kiijv^r^np / d 2 q 



p#0 



(6.32) 



il P \ ( 

q (iS(Op) 



v F q 



x Ira , 



Since the poles in q x from the fermionic Green's func- 
tions are all in the same half-plane, one expects that 
q x yp ~ (7J7) 1 / 3 and thus dominates over the q 2 term, 
which in turn can be neglected in the fermionic Green's 
functions. It then becomes possible to perform the inte- 
gral over q y , which gives 



II A (0,T) = 2iN 



mg 2 T y^ _f2p_ 



n 3 jv F ^-i \Cl„ 



p#0 



Jo 



.1/3 



dq x 



q^og\q x \ 



(iE(Qp) - v F q x ) 3 



(6.33) 



Integrating over the half-space where there is no triple 
pole, we find that the integral is determined by the 
branch cut in log \ q x \. Evaluating the integral we obtain 



dg x ql log \q x \ 
(iY,(tt p ) - v F q x 



- = signup -J log 
|3 Vjr 



E, 



|E(n p ) 



(6.34) 
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FIG. 7: Three-loop contribution to the fermionic self-energy 



Thus 



n A (o,T) = 



2A rag 2 
3 ir 2 ~fVp 



p#0 



Ep 



(6.35) 



To perform the summation over p, we notice that when 
the summand does not depend on p 



A/T 
-A/T 



2AK 



(6.36) 



is independent on T. Then the same sum but without 
p = term will be 2AA — AT. Using this, we obtain with 
logarithmic accuracy 



-Tlog 



T 



(6.37) 



where dots stand for 0(T) terms. Substituting this into 
(|6.35l) . we obtain the final result: 



2kl T , ( E F 

n.(o,T) = ^-fog(^ 



T 
E~f 



(6.38) 



Although small, because of the prefactor in a, this term 
dominates at low temperature over any regular T 2 term. 
The sign of this T log Ep /T term is opposite to the sign of 
a conventional correction to scaling HMM theory, which 
comes from a constant part of the prefactor for the </> 4 
term. In the HMM theory the temperature dependence 
of the spin susceptibility is bT d+z ~ 2 z, b > 0, which in 
d = 2 leads to a linear in T with positive prefactor. The 
negative sign of non-analytic temperature correction in 
(|6.38l) implies that the static spin susceptibility is nega- 
tive right above the QCP. This is another indication that 
the ferromagnetic QCP is unstable. 



D. Self-energy at the three-loop level 

Finally, we show how one can detect the instability 
of a ferromagnetic QCP from an analysis of higher or- 
der self-energy diagrams. This analysis is complimentary 
to the calculations hat we already done in the previous 
subsections. We show that upon inserting the contri- 
butions from the diagrams presented in Fig. [3] into the 
fermionic self-energy, we obtain series of singular correc- 



tions in powers oiu> m i n /u, where (yuj m i n ) 



1/3 



and 



q m in is the scale at which the static susceptibility Xs(q) 
becomes negative, Eq. Q6.29p . 



To illustrate this, we consider one of the three-loop 
diagrams, represented in Fig. In the case of a spin 
interaction, we obtained a finite result after collecting 
the various diagrams, so we restrict ourselves here with 
just one of these contributions. The analytic form of the 
diagram in Fig. writes: 



E 3 (w) ~g J dujid 2 qi 



A(q 1 ,u 1 ) 



- uji) - v F q lx 

\ 2 
1 \ 



2m 




(6.39) 



where A(qi,u>i) is the dynamic fermionic bubble given 
by: 



A(q 1 ,uj 1 ) ~ Ng 2 / dFkdQ. / d 2 q 2 du 2 



1 



1 



i 



7l"2| 
92 



iE(f2) - e fe iE(fi + oJ\ + lo 2 ) - (-k+ qi +q 2 
1 1 



«E(Q - Wi) - ek+ qi iE(fi + lo 2 ) - £fc- 



Q-2 



Approximating A(qi,u)i) by its singular static part 
H^F and substituting into J7J we obtain: 



E 3 (w) ^gy/kp / dw^qx 



qf + y\ui\ 

3/2 

Qi 



iE(w - ui) - v F q lx 



2)ti 



(6.40) 



A simple analysis of this expression shows that the 
dominant contribution to £3(0;) comes from q\ x ~ q\ y ~ 



(70;; 



a/3 



since the integral over q\ x is determined by the 



branch cut in the bosonic propagator. One can then 
safely drop the quadratic term in the fermionic propa- 
gator and integrate over the angle 8 between kp and qi . 
This leads to: 



E 3 (w) 



'0 



dqi 



9/2 



(qf+yuJi) 



g~\fk~F f u du>i 



VFy/l Jo V^l 

Collecting the pre- factors, we find: 

E 3 M ~ (gu) 1/2 ■ 



(6.41) 



(6.42) 



We see that the non-analyticity in the static spin sus- 
ceptibility feeds back into the fermionic self-energy lead- 
ing to a contribution from the three loop self-energy 
whose frequency dependence is more singular than the 
uj 2 / 3 dependence that we obtained assuming that the 
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static susceptibility is regular. Comparing these two con- 
tributions, we see that they become comparable at a fre- 
quency: 



1/3 2/3 



Q 3 

7 



p2 

(6.43) 

9 



where q mm is given by 



Parametrically, q m in is 



not small since q m { n ~ fepj and w m i„ ~ Ep/a is larger 
than However, q TO i n ~ 0.06fc.F is small numerically 
so that uJmin is four orders of magnitude smaller than 
Ep/a. 



E. Summary 

To summarize, we found that the Eliashberg theory for 
an SU(2) symmetric ferromagnetic QCP has to be ex- 
tended to include extra singular terms into both the spin 
susceptibility and the fermionic self-energy. These terms 
originate from the "anti-Migdal" processes in which slow 
bosons are vibrating near the fermionic mass shell. Phys- 
ically, these extra processes originate from the dynamic 
long-range interaction between fermions, which is still 
present at the QCP despite the fact that fermions are no 
longer good quasi-particles. 

We demonstrated that these extra non-analytic terms 
can be understood in the framework of HMM </> 4 theory 
of quantum criticality. We showed that the prefactor for 
the 4> term is non-analytic and depends on the interplay 
between momentum and frequency. The non-analytic 
bosonic self-energy is the feedback from the non-analytic 
<j) A term to the quadratic 4> 2 term. 

We found that these extra terms in the Eliashberg the- 
ory make a ferromagnetic QCP unstable below a certain 
momentum/ energy scale. We detected the instability by 
analyzing the momentum and temperature dependence of 
the spin susceptibility, and also the fermionic self-energy 
at three-loop order. 



VII. STATIC SUSCEPTIBILITY IN THE 
NON-SU(2) SYMMETRIC CASE 

The problem of fermions interacting with bosonic col- 
lective modes with a propagator similar to the one we 
considered, is quite general, and one can wonder to what 
extent our analysis for a ferromagnetic case can be ex- 
tended to the case of an interaction with charge fluctua- 
tions, a nematic QCP, or a ferromagnetic instability with 
Ising symmetry. 

The essential difference between these cases and the 
ferromagnetic one lies in the symmetry of the order 
parameter. In an SU(2) spin-symmetric ferromagnetic 
case, the order parameter (magnetization) is a three- 
dimensional vector, while in the other cases, it is a scalar. 

As we already discussed, the Eliashberg theory, and the 
analysis of its validity at the two-loop level can be car- 
ried out equally for systems with vector and with scalar 



order parameter: the only unessential difference is in the 
numerical prefactors. On the other hand, the evaluation 
of the corrections to the static susceptibility leads to dif- 
ferent results for scalar and vector order parameters, as 
we now demonstrate. 



A. Ising ferromagnet 

Consider first the situation of a magnetically-mediated 
interaction, where we change the spin symmetry of the 
bosons from SU(2) to Ising. The Ising case was argued 
to be relevant for metamagnetic quantum critical points 

M 

The use of Ising spins doesn't change the expression for 
the Green's functions but replaces the Pauli matrix er at 
the fermion-boson vertex by a z . As a consequence, the 
computations performed for the SU(2) case still hold, but 
the interplay between different diagrams changes because 
of a change in the numerical prefactors. In particular, 
instead of Eq. I|6.12f> we now have: 



Ising -plsing ^ 



1 1 — x 2 



r 



Ising -plsing 

J- A 



= o. 



(7.1) 



Under these circumstances, the non-analytic contribu- 
tions from the diagrams in Fig. [5] cancel each other out. 
As a result, the static spin susceptibility remains analytic 
and scales as x~ 1 ((j') ex q 2 with a positive prefactor at the 
QCP. 

This result can be extended to the case of a nematic 
instability, following the same arguments. 



B. Charge channel 

For a charge vertex, one has to replace the Pauli ma- 
trices present at the vertex by Kronecker symbols 5 a p. 
We then have: 

^Charge -pChargc \ X X X X 1 

I 1 = 1 2 d a/30/3 7 1 SOSa = 2, 



-pChargc -pChargc 

1 3 " 1 4 



Substituting these r chargc into the expressions for II, 
we find that the diagrams II ! and H 2 , as well as ±I 3 and 
n 4 cancel each other out. This leaves only regular q 2 
term in the static charge susceptibility. 



C. Physical arguments 

The cancellation of the non-analytic terms in the 
charge susceptibility is not a direct consequence of the 
conservation laws. These laws impose constraints (Ward 
identities) on the behavior of the susceptibilities in the 
opposite limit q — 0, u ^ (Xc{q — 0, u>) vanishes as a 
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uniform perturbation cannot affect a time independent, 
conserved quantity). 

Instead, the absence of the non-analytic terms in the 
charge channel is related to the fact that this suscepti- 
bility measures the response of the system to a change in 
the chemical potential. We showed that the origin of the 
singular behavior of the static spin susceptibility lies in 
the Landau damping term in the bosonic propagator (see 
(|6.24l) . (|6.28l) ). The Landau damping does not depend 
in a singular way on kp (i.e. on the density of electrons), 
and therefore there is no singular response of the system 
to a change in the chemical potential. 

Conversely, the effect of a magnetic field on 
the Landau damping is singular: for a fermionic 
bubble with opposite spin projections of the two 
fcrmions, n± ~ |f2|/-\/0 2 + (vpl) 2 is replaced by 
\fl\/x/(n + 2i f i B H) 2 + (v F l) 2 in the presence of a small 
magnetic field H. As a consequence, if the susceptibility 
is evaluated as a second derivative of the thermodynamic 
potential with respect to H. Taking the second deriva- 
tive of n± over H and setting H = later, one obtains, 
for vpl » f2, a non-analytic d 2 H±/dH 2 ~ Sl/Z 3 . This 
non-analyticit y g ives rise to q 3 ^ 2 term in the static spin 
susceptibility |29L l3?|. For an Ising ferromagnet, this ef- 
fect does not exist as there are no bubbles with opposite 
spin projects in the theory. 

The above reasoning shows that non-analyticity ap- 
pears in the spin response but not in the charge one. To 
further verify this argument, we computed the sublead- 
ing, three-loop diagrams for the charge susceptibility and 
found that the non-analytic contributions from individ- 
ual diagrams all cancel out. We present the calculations 
in Appendix lEl 

The same argument holds for the self-energy at the 
three-loop and higher orders. The singular to 1 / 2 term ob- 
tained in Eq. I|6.42[l appears in individual diagrams, but 
in the case of a QCP in the charge channel (or an Ising 
QCP in the spin channel) the total singular contribution 
cancel out. 



D. Summary 

To summarize, the extra singular additions to the 
Eliashberg theory are specific to SU (2) spin case and all 
cancel out for the charge QCP, the gauge-field problem, 
the nematic QCP, and the spin QCP for a scalar (Ising) 
order parameter. 



VIII. CONCLUSIONS 

We have constructed a fully controllable large- A quan- 
tum critical theory describing the interaction of fcrmions 
with gapless long-wavelength collective bosonic modes. 
Our approach, similar but not identical to the Eliashberg 
theory for the electron-phonon interaction, allows us to 



perform detailed calculations of the fermionic self-energy 
and the vertex corrections at the QCP. 

We constructed a controllable expansion at the QCP 
as follows: we first created a new, non Fermi-liquid "zero- 
order" theory by solving a set of coupled equations for 
the fermionic and the bosonic propagators, neglecting the 
vertex corrections as well as the momentum dependence 
of the fermionic self-energy, and then analyzed the resid- 
ual interaction effects in a perturbative expansion around 
this new zero-order theory. 

We have proved that this approach is justified under 
two conditions: (i) the interaction g should be smaller 
than the fermionic bandwidth (which we assumed for 
simplicity to be of the same order as Ep) 1 and (ii) either 
the band mass raj should be smaller than m = pp/vp, 
or the number of fermionic flavors A should be large. 
When A = 0(1) and ms ~ m, the corrections are 
of order 0(1). We found that the corrections that are 
small in gjEp come from bosons near their resonance, as 
in the Eliashberg theory for the electron-phonon inter- 
action. The corrections small in ms/(Am) come from 
bosons for which one of the momentum component (the 
larger one) is near the bosonic resonance, while the other 
component is close to the fermionic mass-shell. 

For an SU(2)-symmetric quantum critical point to- 
wards ferromagnetic ordering, we found that there ex- 
ists an extra set of singular renormalizations which come 
from bosons with both momentum components vibrating 
near the fermionic mass-shell. These processes can be 
understood as a consequence of an effective long-range 
dynamic interaction between quasi-particles, generated 
by the Landau damping term. These singular renor- 
malizations are not small and have to be included into 
the Eliashberg theory. They give rise to a negative non- 
analytic q 3 / 2 correction to the static spin susceptibility, 
signaling that the ferromagnetic QCP is unstable. 

We also demonstrated that the non-analytic g 3 / 2 term 
can be understood in the framework of the <f> theory of 
quantum-criticality. We showed how the effective long- 
range dynamic interaction between fermions affects the 
structure of the 4 theory, once fermions are integrated 
out: we found that the prefactors of the </> 3 and c/> 4 terms 
appearing in the effective action are non-analytic and de- 
pend on the interplay between the typical external mo- 
mentum and frequency. 

We showed that the non-analytic corrections to the 
bosonic propagator are specific to the SU(2)-symmetric 
case when the order parameter is a vector. For systems 
with a scalar order parameter, like a QCP in the charge 
channel, a nematic QCP, or a ferromagnetic QCP with 
Ising symmetry, the q z l 2 contributions from individual 
diagrams cancel out in the full expression of the suscep- 
tibility. 

The consequences of the instability of the ferromag- 
netic QCP still needs to be fully understood. There 
are two possible scenarios for the behavior of the sys- 
tem: either the ferromagnetic transition becomes first 
order [30l l39| , or the instability occurs at a finite q, lead- 



24 



ing to a second order transition towards an incommensu- 
rate state |4jj The full analysis of these two scenarios is 
clearly called on. 
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APPENDIX A: BOSONIC SELF-ENERGY 

In this section, we compute the bosonic self-energy at 
the one-loop level, in the case of both free and fully renor- 
malized fcrmions. We prove that for an external bosonic 
momentum on the bosonic mass-shell, this self-energy be- 
comes independent on the actual form of the fermionic 
self-energy, and reduces to the usual Landau damping 
term, whereas an extra term has to be included if this 
same external momentum is on the fermionic mass shell. 

After performing the sum over spin matrices, we are 
left with the following expression: 

IT(q,ft)=2iV5 J ^G(k,o;)G(k+q, ll ;+n). (Al) 

Introducing the angle 9 defined by ek+q = £fc + v F qcos9, 
this writes: 

n( q , n) = jvf^ f duo de k do 1 

47T 5 J i(cj + E(a)j) - e k 



1 



i(u> + f2 + S(oj + Q,)) — e k ~ v F q cos ( 



(A2) 



Proceeding with a contour integration over efe, and 
noticing that u> and S(a>) have the same sign, we get: 



n( q ,n) = 



gm 
2tt~' 



+ 00 p2n 

duj I d9 (9(uo + n) - 9(uo)) 

— DO JO 

x \ (A3) 

i(Q + T.(0, + uj)-T,(Lj))-v F qcos9' 

Performing the integration over 9, and rearranging a 
little bit the integration over f2, we are left with: 



mg 



n(q, Q) = N — — / dio Sign(ft) 

1 

X 



(A4) 



We know from direct perturbative calculation that the 
fermionic self-energy goes like = uj^ 3 uj 2 / 3 , where 

luq ~ ag. It follows that if the external bosonic mo- 
mentum is on the bosonic mass-shell, i.e. if q ~ fi 1 / 3 , 



it dominates over the frequency-dependent term in the 
integral, so that: 



tt/ ^ ,r m 9 f Q , Sign(fi) 
n(q, O) = N — - / du & v ; 

7f Jo v Fq 



N 



mg |0| 
irv F q 



(A5) 



We recover here the expression of the Landau damp- 
ing term with a prefactor depending on the details of 
the interaction considered in our model. This result is 
independent on the fermionic self-energy provided that 
v F q^> (ft + £(£))). 

However, if the external bosonic momentum is on the 
fermionic mass shell, i.e. if v F q ~ £(f2), the frequency- 
dependent term can no longer be neglected. Defining 



one then has: 



n(q,o) = A^ 

TT 



x / dz 
'o 



V Xm) 2 + ( zV3 + (! - z ) 2/3 ) : 



(A6) 

This formula is of limited use as it is modified by vertex 
corrections. For the calculations of the static spin sus- 
ceptibility, we will actually only need the leading 0(1/ q) 
and the subleading, 0(l/q 3 ) terms in H(q ) Q) We will 
show in Appendix C that these two terms still can be 
evaluated without vertex corrections (see Eq. (|C20|) be- 
low). Expanding (|A6(I in 1/q to first two orders, we find, 
at < u>o, 



nq,fl=7 U i A 2 

q \ v'q' 1 



(A7) 



where c ~ 1.19878. To simplify writing, we will be us- 
ing the expression for the polarization operator with f2 
dependence in the denominator, i.e., plug E 2 term back 
into denominator, i.e., 



n(q,Q) = iV^ 



\n\ 



^ ^(v F qf + c 2 ^(nf 



(A8) 



but one should keep in mind that we will actually only 
need 1/q and 1/q 3 terms in the calculations. 



APPENDIX B: FERMIONIC SELF-ENERGY 

In this section, we compute the fermionic self-energy 
at the one, two and three loop levels, for an external mo- 
mentum taken to be on the Fermi surface. We also ana- 
lyze the momentum dependence of the one-loop fermionic 
self-energy. 
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1. One-loop 

a. At the Fermi level 

After summing over the spin matrices, the fermionic 
self-energy at the one-loop level is given by: 

E(w) = Zig 2 J ^^G(k P + q,a; + n)x(q,fi) 
3z.g 2 



(2tt) 3 



dn dq d9 



Xo 

e - 2 + g2+7 M 



q 



i{u> + VL + E(w + fi)) — VFqcosO 



(Bl) 



where we defined as the angle between k and q, and 
considered an external fermionic momentum k ~ fc^?. 
The i prefactor comes from the convention G -1 = Gq 1 + 
iS. 

Since the pole in g from the fermionic propagator is 
in a definite half-plane, the integral in q is dominated by 
poles coming from the bosonic Green's function, so that 
one can perform the integral over the angular variable 
and simplify the result as follows: 



E(w) 



3g 2 



dfl dq 



Ao 



q 2 Sign(w + fi) 
V(uf9) 2 + + n + E(w + «)) 2 



(2tt) 2 







1 



(B2) 



Defining the new variables z — — and u — 2 s2 , this 
leads to: 



2 /•! 



E(w) 



/ (iw 



3ff _w 

27r 2 7 i>f, Jq Jo z + a« + 6m 3 



> (B3) 



1 uo 2 
where we have used a = - and b = — ~- 

Let's denote by X the last double integral, and define 



the new variables y = J -^z — (luj(, )z and t 



&-u. Now I reduces to: 



— \ 2 r" e d f +o ° dt 1 

W Jo V Jo 1 t 3 +t + y 



Dp 



(B4) 



Substituting this back into our expression for the self- 
energy: 



E(w) 



27rm£ : 



(B5) 



7r.r 
2 



with the following asymptotic behavior: h(x — > 0) 
and /i(a; — > oo) = . 

In the two regimes we are interested in, this last result 
can be rewritten as: 

v , v J Aw for > 1, , . 

S M = \ 1/3 2/3 f c -i n ( B6 ) 
I w a; ' for £ — > 0, 

where the constant prefactors are defined as: 



A = Ali 



UJ 



47T Vf ' 

3^3 f 
8ir 3 jv F ' 



(B7) 



6. Momentum dependence 

For definiteness, we set = 0. We first compute 
the momentum-dependent part of the one-loop fermionic 
self-energy at zero external frequency: 

t d 2 a dO 

E(k,w = 0) = iig 2 J 7 ^ r G(k+q,%(q,U) 
3ig f dnd 2 q q 



(2^) 3 7 iE(fi) - e fc+9 7l«l + 9 3 

Expanding efe +g = + vpq cos 9 and integrating over 
we obtain 



E(k,0) 



4tt 2 



dfi Sign(fi) 



q 2 dq 



Q 3 + lM 



(B9) 



((^g) 2 + (E(r!)+ie fe )- 



1/2' 



At ej, = 0, the integral vanishes by parity. Expanding to 
linear order in efc we obtain: 



(BIO) 



E(k, 0) = -Me k x / dft E(fi) x 

<f<k 

(g 3 + ( 7 |fi|))((« F q) 2 + E(0) 2 ) 3/2 

Simple estimates show that the result depends on 
the interplay between (■yfl) 1 ^ 3 and E(f2)/«ir. Introduc- 
ing the scale WMax, defined as the frequency at which 
(7fi) 1 / 3 = S(fl)yVF (see (|5.3|) . and rescaling variables as 
<7 = (7WMax) 1 ^ 3 y, H = WMra^, we re-write (|B 10|> as 



L(k, 0) = — 3i£fc x „ Q — j 



2tt 2 w 3 



«f7 



(Bll) 



where 



1=1 dx I dy 



o Jo (x 2 + y 2 ) 3 / 2 (x + y 3 ) 



1.3308 



(B12) 
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Substituting WMax = (t^f) 1 ^ 2 , we obtain l|5.27[l . 

To obtain the frequency dependence of the fermionic 
density of states at small w, , we have to evaluate the sec- 
ond order term in e k at the mass shell, where e k = iE(cj) 
and convert the result to real frequencies. We therefore 
will keep both to and e k finite, and use 



(2tt) 2 



1 



lM + q 3 i±{uj + Q) -€ k+q 
(B13) 

Writing, as before, e k + q = £fc + «f9 cos # and integrating 
over we obtain 



3g f dilq 2 dq 



sign (a; + Q) 



We assume and then verify that the internal VFq are still 
larger than E(w + 0) and e kl and expand 



E(k,w) 



3.9 



8tt 2 4 



<if2sign(u> + £1) 



tig 



g( 7 |n|+9 3 ) 



T,(uj + n)+ie k ) (B15) 



The lower limit of the momentum integral is (E(w + f2) + 
ie k )/vF- At the mass shell, ie/c = E(w). Substituting, we 
find 



E(k,w) 



3.9 



8tt 2 4 



dSl sign(o; + Q) 



9( 7 |«|+9 3 ) 



E(w + ft) - E(w) (B16) 



We will need the contributution which is confined to 
ft ~ to. The contributions from \Q\ 3> |w| diverge in our 
expansion procedure, and account for the regular 0(e k ) 
and 0(lo) terms in the self-energy. The last term is even 
smaller in a than the regular 0(e k ) term and is totally 

irrelevant. As ft ~ u is small, E(ft) w E(ft) = ft 2 / 3 ^,!/ 3 . 

Because of the sign factor in the numerator of l|B14|l . 
there are two distinct contributions from Q ~ w. For 
both of them, the momentum integral is logarithnic (this 
justifies the expansion) and yields (1/3) log(wi/cj), where 
u>\ ~ N 2 E F /g. The first contribution comes from |ft| < 
| a; | and to logarithmical accuracy is 



87T 2 «|,7 



M dft 

sign(w) / T7T7 



ft 



E(k,w), 



Rescaling the frequency, we obtain from l|B17(l 



x(£(M + ft)-E(M)) 2 logfi(B17) 



E(k,u/U 



4-7r 2 'i; 3 - 



V- E(w)|E(w)|log-i. (B18) 



«f7 



where 



/i = i/ 1 i ^((l+-) 3/2 -l) 2 = 0.254 (B19) 



Another comes from 101 > \u)\, and is 



„,, , <7S(cj) , wi Z" 00 <ift 



4tt 2 4 7 b | W |7 H 
x (E(|w| + ft) + E(ft - - 2E(ft))(B20) 



Rescaling, we obtain 



-< k -\>" : i £ Ml s ^)l lo gn- (B2i) 



47r 2 w|,7 



where 



Z 2 = / ^ ((1 + x) 2 / 3 + (z - l) 2 / 3 - 2a; 2 / 3 ) 
= -0.195 (B22) 

Combining E(k, u>)a and S(k, we obtain 
0.455 



E(k,w)= 



47T 2 w|,7 



E(w)|E(w)|log^. (B23) 



Substituting the result for 7 in this last expression, we 
obtain (IQTl). 



2. Two-loop 

We compute here one of the contributions to the two- 
loop self-energy, given by Fig. ^ This contribution orig- 
inates from the insertion of the vertex correction into the 
Eliashberg self-energy. 

We have: 

E 2 (w) ~ g 2 J dwid 2 <ji J duj 2 d 2 q 2 x(<li,^i)x(<l2, u 2 ) 

xG(k F + qi,w + LUi)G(k F + q 2 , u + uo 2 ) 
xG(k F +qi + q 2 ,uj + loi +lu 2 ), 

(B24) 

which gives, once we replace each propagator by its full 
expression: 



E 2 (w) ~ g / du)\d qi / du 2 d q 2 



J2„ 91 



92 



1 



7|wi| +ql 7|w 2 | + 5 2 
1 



iE(w + wi) - «f9i* it{io + u 2 ) - v F q 2x - ^ 
1 

iE(w + Wi + W 2 ) - UF9lz - - - ,777- 



(B25) 



where we use the shorter notation E(w) = a; + E(u>). 
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Integrating successively over q\ x and q 2x , closing each 
contour on the upper half-plane, one has: 



S 2 (w) 



m B g 



duoi I du>2 uj\, w 2 ) 



dqiy J dq 2y 



Wy\ 



\<liy\ 3 +7|wi| \q 2y \ 3 + j\uj 2 \ 



(B26) 



where £ = tub 



£(w + cji) + + W2) — S(w + uii + W2) 



(7 2 wiw 2 ) 



and Q(cj,u)i,u>2) comes from the choice of a contour for 
the integration and is given in our case by: 

®(w,u)i, uj 2 ) = {6{u) + ui) — 9(uj + uji + lo 2 )) 

x (e(u + uj 2 ) - 0(S(w + LUX+ lu 2 ) - E(w + wi))J . 

(B27) 

It is convenient at this stage to rescale the perpendic- 
ular components of the bosonic momenta. Introducing 
x = 9iy/(7ki|) 1/3 and V = q2 V /(l\u 2 \) 1/3 , we obtain: 



£ 2 (co>) 



dcu 2 / du>i / dxdy 



6(w,wi, L0 2 ) 



o (7 2 wiw 2 ) 2/3 
xy 



x 2 y 2 + ( 2 (1 + a; 3 ) (1 + y 3 ) ' 

(B28) 

where we rearranged the double integral over x and y to 
make it real. 

Since all internal frequencies typically go like u, the 
typical value of £ is given by the small parameter given 
in (15.13(1 . Expanding the double integral for small values 
of C, the leading contribution from the integral over x 
and y reads: 



1 



xy 



„ ^ x 2 y 2 + ( 2 (l + a;3) (1 + 2/3) 



log 2 C. (B29) 



If one considers now free fermions, it becomes possible 
to reduce the expression of the two-loop self-energy to: 



™ 2 « 2 /■<" log 2 ( , 2 mB ^2/ 3 

S free (w) „ 1 dul2 r duli _ VW^WS 



V F JO Ju>-uj 2 

— ^5— u> / dz 2 / dzi 

IVp Jo Jl-z 2 

(3 2 u; log 2 lo, 



^ A LO\U} 2 



l 0g 2 ( ™>) 



ZlZ 2 



(B30) 



where we only kept the leading contribution in the last 
expression, and [3 — ms /mN is one of the small param- 
eters defined in the text. 



Now, in the case of dressed fermions, we need to take 
= luq^uj 2 / 3 . The procedure is identical to the free 
fermion case, but the final expression is a bit more com- 
plicated: 



■r-n 2 n 2 , P L Z" 1 I 

.Dressed/, A ^ 171 9 1/3 2/3 ' ■ ' 1 



? 2 ^0 



1 f l 

dz 2 I dz\ - 

Ji- Z2 z x z 2 



X (l-Z 1 ) 2 / 3 + (l-Z 2 ) 2 /3 + ( Zl+Z2 -l) 



2/3 



^2/3 



x loe 



(1 - Zi) 2 / 3 + (1 - z 2 ) 2/3 + ( Zl +z 2 - l) 2 / 3 



7 ^/3 1/3 1/3 
1/3 z l z 2 



(B31) 

Expanding the log^, one is left with a double integral 
that only contributes as a numerical prefactor, and the 
dominant term is then given by: 



E 2 (w)~Ei(w) /? 2 log 2 /3, 



(B32) 



where Y>i(uj) = lUq^uj 2 / 3 is the self-energy in the Eliash- 
berg theory. 



3. Three loop 

We now turn to the computation of the three-loop self 
energy. We are only interested here in one diagram, given 
in fig. U| where we try to analyze the feedback of the non- 
analytic susceptibility into the the higher-order diagrams 
for the fermionic self-energy. For spin interaction, there is 
no cancellation between different diagrams for the static 
susceptibility, which justifies that we restrict ourselves to 
just one contribution. 

The analytic expression for this diagram is: 



E 3 (w) ~g J duid 2 qi 



A{qi, uji) 



i£(w - uji) - v F q lx - 



1 



2 , TkjJ / ' 



(B33) 



where A(qi,LUi) is the factor from the fermionic bubble: 



A(q 1 ,uj 1 ) ~ Ng 2 I d 2 kdQ, \ d 2 q 2 dbj 2 

J J a 2 + - 



1 



12 
1 



7l"2| 

92 



i£(Q)-e fe iE(ft + wi + lo 2 ) - e k+qi+q2 
1 1 



iE(f2 - uji) - e k+qi iE(fi + lo 2 ) - e k+q2 ' 
Approximating A{qi,oji) by its singular static part 

3 /2 1 

Qi ykp and substituting into the expression of E3 we 
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obtain: 



£3(0;) ^g\fkp J duj\d' 



'7i 



3i 
3/2 



iE(w - wi) - ^gi^ 



2 m B 



(B34) 



The integral over gi^ is determined by the branch-cut in 
the bosonic propagator and one then expects that this 
very integral is dominated by q\ x ~ (7^1 ) 1/>a . It follows 
that the term in q\ x dominates inside the fermionic prop- 
agator allowing us to neglect the curvature term. Defin- 
ing the angle 9 between k F and qi , and integrating over 
it, this leads to: 



£ 3 (w) 



kp f dujidqi sign(u; — u>i) 



9/2 
9i 



s("-"i) 2 (gf + tI^iI) 2 



(B35) 

Since the dominant contribution comes from q\ ~ 
(7( J Ji) 1 / 3 , one can neglect the fermionic self-energy in the 
denominator. This in turn allows to simplify the fre- 
quency integral, which then writes: 



S 3 (w) 



VF Jo 



duji J d 9ij^ 
du>i 



7/2 



) 



(B36) 



where we introduced z = q\/ (7W1) 1 / 3 , so that the integral 
over z just contribute to the numerical prefactor. 
Collecting prefactors, one finally has: 



S 3 (w) 



gu. 



(B37) 



APPENDIX C: VERTEX CORRECTIONS 

In this section, we compute the various vertex correc- 
tions analyzed in the text. 



1. q = Q = 

Consider first the simplest 3-leg vertex, with strictly 
zero incoming frequency f2 and momenta q, as presented 
in fig. Its analytic expression writes: 



Ag 
9 



q=Vl=0 



J dujd 2 p G(k F ,uj) 2 x(p,u;) 
dtod 2 p 1 



7H 



T+P 2 (iE( 



u) - v FPx - ^ 



where we defined £(w) = lo + E(oj) and we have chosen 
hp along the x axis. 



Since both poles coming from the fermionic Green's 
functions are in the same half plane, the integral over q x 
is finite only because of the branch cut in the bosonic 
propagator. Since at the branch cut p x and p y are of 
the same order, we can drop the curvature term in the 
fermionic propagators and introduce polar coordinates 
for the internal bosonic momentum. Defining 9 as the 
angle between k F and p, and integrating over it, one 
has: 



Ag 
9 



9 



q=Q=0 



dujdp 



p 2 |EH| 



p 3 + 



3/2 • 



( VF p) 2 + EH 2 j 

( C1 ) 

Introducing the frequency LJ max up to which bosons 
are slow modes compared to fermions, i.e. up to which 
(7dj) 1//3 > T,(lo)/v f , one can split the frequency integral 
into two parts, and define in each of them the reduced 
momentum z = p/Min((7w) 1 / 3 , T,(lo)/v f ) so that the in- 
tegral over z only contributes to the numerical prefactor, 
leading to: 



Ag 
9 



EH 
ivpuj 



9 fv 

7<4 



(C2) 



One can easily make sure that the frequency w max at 
which (7a;) 1 / 3 = T,(lu)/v f exceeds ojq, such that one 
should use E(w) = lo to find Lo max . Substituting, we 
obtain u max ~ (NgEp) 1 / 2 , and 



9 



fa 



(C3) 



q = Q = Q 



2. q = 0, n finite 

Considering the same vertex, now with a finite external 
frequency, one has: 



Ag 
9 



9=0 



~ g 2 dujd 2 p G(k F + p, lo + Q) 
a J 



G(k F +p,w)x(p,w) 
dcod 2 p 1 



P 2 iV{u)- v FPx -^ 
1 



, (C4) 



it(u + n) - ufpx - ^ 

where we chose the x axis along k F . 

From the pole structure in p x of this expression, one 
expects two contributions to this integral. A first term 
comes from the branch cut in the bosonic propagator, 
however this contribution ultimately leads to the same 
result as the q = Cl = vertex, up to small corrections 
from the finitcncss of Q. The second contribution arises 
from taking the poles in the fermionic propagators. At 
zero external frequency, these two poles were in the same 
half- plane of p x , so we could close the integration contour 
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over a different half-plane and only consider the contri- 
bution from the branch but in the bosonic propagator. 
At a finite fl, there is a range where u> and cu + tt have dif- 
ferent signs, and the two poles are in different half-planes 
of p x . The result after integration reads: 



Ag 
9 



q=o,n 



9 f n f 
— cLo dp v 

VF Jo J 



\Py\ 



iH + \Py? 



X — 



t(si-w)+£(wy 



(C5) 



where we slightly rearranged the frequency integral. 
Performing the integration over p y , we are left with: 



Ag 



<?=o, n 



7 1/3 «f Jo 



duj — 



,-1/3 



E(Q - u) + E(w) 



(wo74) 1/3 
Const., 



(C6) 



where we assumed that il is small, i.e. E(f2) = Wq^ 3 !} 2 / 3 . 
This vertex thus reduces to a numerical constant, that 
does not contain any small parameter. 

3. q finite, f2 = 

Conversely, the same vertex taken at finite external 
momentum q, but zero external frequency writes: 



Ag 
9 



q,n=a 



g J dud p G(k F + p + q, lj) 

G(k F +p,w)x(p,w) 
dujd 2 p 1 



^+P 2 iZ(u)-v F p x 



P 1 17 "rrx 2m B 

1 



iE(w) - VFQx - VFPx 



2rriE 



IvPv 
m B 



where p x is the projection of p along k F . 

Like its q = counterpart, this vertex is characterized 
by poles in p x from the fermionic propagators lying in 
the same half-plane. The only nonzero contribution then 
comes from the branch cut in the bosonic propagator. 
At the branch cut, p x <~ p y which allows us to neglect 
the quadratic curvature terms in the fermionic Green's 
functions. 

This makes possible a direct integration over p y . This 
integral can be separated from the rest of the expression, 
and reads: 



dp. 



P 



7|w| +p 3 



\Px\ 

7M Jo 



dz 

-co 

+ 7T/2 



VTTz 1 



( 1 + z 2)3/2 +7 JM 



Ip*I 3 



du 



^ + (cosu) 3 ' 



(C7) 



where we successively defined z — p y /\p x \ and z = tanu. 

This last integral can be approximated by its asymp- 
totic form, namely: 



/ 



dp y P 

7M +p s 



log 



Px 

3 



1^, if |p x | 3 > 7 w 

4tt 1 _ 1 Pi 
3 V / 3(7|"I) 1/3 27R 

+ (log2-i)^, if b,| 3 «7^ 

(C8) 

The only non-vanishing contribution once we take the 
integral over p x comes from the log term. After expand- 
ing in q x , this contribution reads: 



Ag 
9 



q,n=o 



Ag 
9 



gvF 



q =n=o 

lv\l dPx 



Pi iogPx 



(it(uj) - v F p x ^j 



f (C9) 



Defining the scaled momentum z = vfPx/^(^) and 
performing the integration over z, one obtains two terms, 
the dominant one being: 



Ag 

9 



q,Q=0 



Ag 
9 



■ 9 f du> ~ 

I — Tlx / — log «E (u , 

7«F J\vpq m \ W 



q =n=o 

(CIO) 

where the frequency integral runs over |E(w)| > |«f9x| 
since as we expanded in VFqx, we assumed that it was 
smaller than |E(w)|. 

Performing the remaining integral, one finds: 



Ag 

9 



q,n=o 



Ag 
9 



q x 
kF 



l0g|?x|- 



(Cll) 



4. q, Q finite 

Finally, we consider the general vertex where both ex- 
ternal bosonic momentum and frequency are non-zero. 
In analytic form, this writes: 



Ag 

9 



q,Q 



g 2 j dud 2 p G(k F + p + q, u + 0) 

G(k F +p,w)x(p,w) 
I" dujd 2 p l 

9 J lH +p 2 



1 



Py 

2m B 



it{u + ^)-V F q x -V F p x -^-^ 

where is defined as = p • k F . 

Integrating over p x first, there are two contributions. 
One comes from the branch cut in the bosonic propaga- 
tor, and gives similar results to the (q = 0, = 0) and 
the (q finite, O = 0) vertices up to small correction from 
the finiteness of the external frequency. We neglect it 
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here and focus on the other contribution which comes 
from the poles in the fcrmionic propagators: 



A5 

9 





9 f U f 
•> i — / dui I 




Vf Jo J 



\Py\ 



j\Uj\ + \Py\ 3 



-u) + *E(w) - v F q x - ■ 

where the simplification of the frequency integral comes 
from the poles in p x . 

This vertex correction strongly depends on the inter- 
play between the external q x , q y and f2, and is in par- 
ticular quite sensitive to the momentum anisotropy. We 
now analyze the various possibilities. 

For the generic case where q x ~ q v (we use the nota- 
tion q to designate them), one can neglect the quadratic 
term in the fermionic dispersion, allowing to perform the 
integration over p y , leaving us with: 



dto or 1 / 3 



q n VfT' 6 Jo iS(fl - u) + iS(u}) - v F q 

(C13) 

Restricting ourselves to the quantum-critical regime (i.e. 
< ^>o) f° r which E(a>) = uJq^lu 2 / 3 , one has: 



A3 
9 

1 dz 



q,fl 



V F q 



(C14) 



where F{x) = f Q -4^-^ 
asymptotic behavior: 



s(n)_ 

(i_ 2 )2/3 +z 2/3 +m ; has the following 



T{x<^\) = 0(1) 
T{x»l) = O(i) 



(C15) 



If the typical q is on the bosonic mass shell, then q 
^fi) 1 / 3 , and one has: 



AS 
9 



q,Q 



E(Q) 



1/3 



(C16) 



However, we encountered in previous computations 
(e.g. self-energies) that a strong anisotropy can be ob- 
served between the components of the bosonic momen- 
tum, with q y 3> q x . In this case, the full expression of 
the vertex correction is a bit complicated and we choose 
to present here the most relevant case for which the cur- 
vature term dominates over v F q x in the fermionic prop- 
agator. The vertex correction then no longer depends on 
q x and writes: 



A.9 



9 



n />oo 
did I dpy 

q,Q V F Jo JO 



Py 



7cj + pi 



£(fi-w) + E(w) 



X (e(o-,) + eh) 2 + (» 



(C17) 



Defining u = py/^ui) 1 / 3 and z = uj/fl, it is possible to 
rewrite the vertex correction in this regime as: 



Ag 



G[0 



q,n 



,(70) 



1/3 



Qv 



(C18) 



1 (C12) where /3 = and G(x) is the following double integral: 



9(x) 



du 



1 dz Z- 1 / 3 ((l-z) 2 / 3 + 2 2 / 3 ) 



1 + JO ((l~z) 2 / 3 + Z 2 / 3 ) + 

a; 2 log 2 a; if a; -C 1. 



u 2 z 2/3 



(C19) 



Finally, for the computations of the full dynamic po- 
larization bubble will also need the vertex averaged over 
the directions of q. The generic structure of this vertex, 
which we define as < Ag/g > is the same as in (|C14ll . 
i.e., 



Ay, 

9 ' 



T 



VFq 
E(O) 



F(0) = O(l), T{x^>l) = 



(C20) 



However, it will be essential for our further analysis that 
the expansion of !F(x) at large x holds in odd powers of 
1/x, i.e., T(x ^> 1) = ai/x + a 3 /x 3 + .... In particular, 
there is no term 0(l/a; 2 ), which we found in the polar- 
ization operator without vertex corrections (see 1A7J0 . 



5. 4- leg vertex 

In this paragraph, we compute the renormalized 4-leg 
vertex T 2 (q, presented in Fig [3 which reads: 



r a(g ,n) ~ f J *> J d 2 P Xs (°±^,E+S 

Q — uj q — p\ /f2 + w P + q 



2 ' 2 



G 



,k F + 



(C21) 



xG ( — — ,k F + — 



Performing the integration over p x , projection of p 
along kp, we obtain: 



2 r n 



5_ 

w F 7 7(fl - + (1 2 +p 2 y - 2q v p y ) 3/2 



yj (q 2 + p$) 2 - 4qfo 



' 7 (fi + w) + ( 9 2 + p 2 + 2q y p v YI* ' 



(C22) 
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It is convenient at this stage to define the reduced vari- 
ables z — to/fl and y = p y /\q y \: 



( 7 n) 



1/3 



dz dy 



1 



i[(l-z) 2 /3 + (l + z )2/3] _ 



^(i + z) + ii+ y r 



s(n) 



l+jr 

/9 



2/3 



(C23) 



where we assumed that < c^o, and we used that in all 
our computations, the perpendicular component of the 
bosonic momentum is always either dominant or compa- 
rable to the parallel one, so that one has q y ~ q. 

Comparing this renormalized vertex with the bare one 



given by Tx(q, Q.) 



i+ 



two: 



Tin I i 

\ly \ 3 



r 2 ( 7 ^) 1/3 



one has for the ratio of the 



dz dy 



o t$t(1-*) + \l-V? 



The remaining two diagrams (their total contribution is 
Hsiq, 0)) cannot be computed explicitly at QCP, and we 
compute them in an approximate scheme. 

In explicit form, the first two diagrams in Fig. [5] are 

given by: 



ni„(g,o) = r a 



(2tt) 6 



d 2 Kduid 2 ldn G(lo, kf 



x G{u + n,k + l)G(uj,k + q)xs(l,ty (Dl) 



U lb {q,0) =r 



9 



d 2 Kdcud 2 ldn G(uj, k) 



(2tt)6 

x G(u> + CI, k + l)G(u) + Cl,k + q + l) 
xG(w,k + q) Xs (l,n) (D2) 

where r a (, are numerical prefactors coming from spin 
summation. Note that for symmetry reasons, one has 
to count the first diagram twice, so that the total contri- 
bution reads: 



U A (q,0) = 2n lo (g,0) + II 16 (<7,0). 



(D3) 



' i [(1 - z)2/3 + (1 + ,)V3] - ^ - 

|i-y 2 l 



l£pr(l + *) + |l + tfl 



3' 
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which only depends on two parameters: the ratios f^y 

\ly\ 3 



and 



In the generic case of an external bosonic momentum 
on the mass shell, i.e. q x ~ q y ~ (7s!) 1 / 3 one has: 



1/3 



(C25) 



On the contrary, for a typical g x ~ S(fi)/uj? and q y 

( 7 fi) 1 /3; 



£2 

Ti 

up to logarithmic factors. 



(C26) 



a. first 



To prove our point, we try to expand the products of 
fermionic Green's functions into a simpler form: 

G(u, k) 2 G(uj + n,k + 1)G{uj, k + q) = 

G(uj,k) 2 G(uj,k + q) G(uj,k)G(uj,k + q) 



G(u + n,k + l)G(u,k + q) 

+ K (i,uj,n) 2 ' 



(D4) 



where k(1, u>, ft)i(S(w + f2) — S(cj)) — Vpl x 

The interest of such a splitting up is that one can re- 
duce this drastically by performing the integration over 
k. In fact: 



de k G(to,k)G(uj,k + q) = 0, 



(D5) 



APPENDIX D: STATIC SPIN SUSCEPTIBILITY 

1. Diagrams 

In this Appendix we present the details of our calcula- 
tions of the singular terms in the static spin susceptibil- 
ity. We discuss in great detail the calculation of the first 
two diagrams in Fig. [5] (vertex and self-energy correction 
diagrams) . These two diagrams can be computed explic- 
itly both away from QCP and at QCP. We labeled the 
total contribution from these two diagrams as 11^(5,0). 



since all the poles in Sk are in the same half-plane. 
For this reason, we are left with: 



n la ( g ,o) = 2r a 



9 I rl 2 1 , , 2 

(2tt) 6 



d 2 Kdujd 2 ldn Xs(l,ty 



G(uj + n,k + l)G(uj,k + q) 



Let's keep this expression as it is for the moment and 
move on to the second diagram. 
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b. Second diagram 

Following the same path, we can rewrite the product 
of fermionic Green's functions as: 

G(uj, k)G(u + fl,k + 1)G(uj + fl,k + l + q)G(u, k + q) = 
G(u, k)G{u), k + q)- G{u, k)G(u; + Q,k + l + q) 



+ 



k(1, u>, O) 2 
G(uj + Sl,k + l)G{u + tt,k + l + q) 

k(1, u, O) 2 
G(u> + fl,k + l)G{w, k + q) 

n(i,to,n) 2 ' 



(D7) 



with the expression of k defined above. 

Once again, the integration over e k may give zero if 
the poles are in the same half-plane, which reduces our 
previous expression to: 



ni6(g,o) = -r. 



(2tt) 6 

G(u + fl,k + l)G(u, k + q) 

n(i,uj,n) 2 

G{u,k)G(uj + Q,k + l + q) 



+ 



(D8) 



Changing k into k — q in the second part of the integral, 
we have: 



n 16 ( g ,o) = -r 



■9 / j2 tsa, .2 



d 2 Kdud 2 ldn Xs{l,V) 



(2tt) 6 

G(u + Q,k + l)G(u, k + q) 
k(1, u>, fl) 2 
G{u,k- q)G{uj + n,k + l) 



+ 



k(1,u>, f2) 2 



(D9) 



One can then notice that £k- q — £k — vpq cos 9 changes 
to Ek+q if one changes 9 into 9 — ir. This finally leads to: 



n lb (g,0) = -2T b7 ^-^ I d 2 Kdu;d 2 ldn Xs{l,ty 

(D10) 



'(2tt) 6 

G(w + n,k + l)G(u},k + q) 



k(1,U!, il) 2 

^From what precedes, we have: 

n A (g,o) = 2fi-^ni»(9.°)- ( Dn ) 



Spin summation prefactors can be easily computed, and 
are given by: 

| T a = E a ,/3, 7 ,{ <J f/3 (T /37 ' ^fS^Sa = 6 
1 T b = E a ,/3, 7 ,J <J f / 3 <T /37 cr 7'5' T ' 5ct = ~ 2 

This finally leads to: 

n A (g,0) = |ni„( g ,0). (D12) 



2. Away from the QCP 

In the Fcrmi-liquid regime we have 



e fe ) 2 



1 



1 



i(l + A)(w + Q) - e k+ i i(l + X)lu - e k+q 

X 7TT77— W (D13) 

where we used for the fermionic self-energy £(o>) = Xw, 
since we are deep in the Fermi liquid phase in this case. 

Defining cos 9 = and cos 9' = f^p^y , and integrat- 
ing over k and u, one has: 



f 

Jo 



dl I 



n 



£-2 + p + lM + A)n - M cos6») 

1 



i(l + X)Q — vpqcos9' — vfIcos9' 
The integral over 9' then gives: 

n A ( g ,o) = -^^ C do r dn r di 

71" JO JO JO 

Z 

£- 2 + I 2 + 5p ((1 + A)fi + ™ F Z cos 6») 2 

i 



(D14) 



-y/ (vFg) 2 + + A) + w F i cos 0) 2 



(D15) 



It is convenient to rescale the variables at this stage, 
introducing f2' = ( 1+A )° and V = -, so that the previous 

° vpq q 1 

expression reduces to: 



iTtM) 



ANmg 



z,2 



7T 4 U F (1 + A) 2 



nTT f-\-00 pOO 

q\ d9 dn' / dl' 
Jo Jo Jo 



1 

x . 

^1 + (W + il' cos 6) 



(D16) 



where we kept only the leading order in q. 

Defining z and (ft as z cos (ft = V and z sin (ft = Q', one 
is left with: 



+ -yg_ tan ^ g . n ^ + . ^ ^^2 

x 1 (D17) 

W 1 + z 2 (sin + z cos (ft cos 9) 
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Subtracting the constant part 11^(0, 0) (and neglecting 
it), and integrating over z, this leads to: 



ANmg 2 



7r 4 « F (l + A) 2M y £- 2 + ^tan 
1 



tt/2 



cos m sin ( 



x / d6- 



o (sin 4> + i cos (j) cos ( 



(D18) 



The angular integration over 9 can be done explicitly 
and gives: 



L 



d0 



o (sin 4> + i cos <fi cos 9) 2 



sin 0(5 sin 2 0-3). (D19) 



Substituting this into the expression of SHa, we are 
left with the following final result: 

'"'"•'-^"''(s?)' (D20) 

where 7i is defined as: 

■u( \ ^ cos sin 2 0(3 -5 sin 2 0) 

Jq idii cp -\- x 

In the two limits, 

W(0) = i W( a; >l) Rj A (D22) 

As one approaches the QCP, £ gets bigger, and one can 
take the asymptotic form of H(x) for small x: Tt(0) = 
1/3. Rearranging the prefactor for this limit, we are left 
with: 



3. At criticality 

At the QCP, we have: 

U A (q,0) = ^£ I d 2 Kduod 2 l<m ' 



(D23) 



(27T) e 



1 



(i£(w) - e fc ) 2 
1 



i£(w + fi) - e fc+; iS(w) - e fc+9 
1 



/ 2 + 7 ¥' 
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where we considered that, close to criticality, the self- 
energy dominates completely the bare uj term in the 
fermionic propagators. 

i 2 

We expand both energies as Ck+i = ek + vfIx + 
and Ck+q — £fe + vf<7 cos 0, and perform the integration 



over Ek, leading to: 



IQNma 2 [ 2lx f +0 ° f +0 ° f n 
IL A (q,Q) =i—^—^f- I d9 j dl x I dQ da; 

J -oo JO 



(2^) 5 Jo 



dly I 



— OO 

1 



uj) + iE(w) - v F l x - ^ 
1 



— lu) + iTi(u>) + vf<1 cos 9 — Vpl x — 



2mi 
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The integration over l x brings two contributions. One 
comes from the poles in the fermionic propagator, and 
can be neglected here since both poles are in the same 
half-plane. The other contribution comes from the 
branch cut in the bosonic propagator, and since at the 
branch cut q x ~ q y , one can safely drop the quadratic 
term in the fermionic propagators. This allows us to 
integrate over l y . Out of the terms arising from this in- 
tegral, the only non-vanishing ones come from the non- 
analyticities of the integrated bosonic propagator defined 
as: 



dl y x{\ fi) = / dl y 



1 



l 2 + 



-, (D26) 



% /(d f () 2 +c 2 S(^) 2 



We use here the full form of the polarization operator, 
Eq. (|A8|) as we will see that typical vfI x > S(^), and 
typical l y are only larger in logarithmic sense. 

This integral was performed in a slightly different form 
in l|C8|l . but the method is the same: introducing u such 



that tan u — 



Vfly 



f, one has: 



di yX (i,n) 



^(v F l x ) 2 +Y,(Q,) 2 

(v F l x ) 2 + c 2 X(tl) 2 



7^1 

it/2 



du 



cos 3 u — 8 cos 2 u + e 



, (D27) 



where we introduced S — c s ^ ^" Fil jin, c s ^ - and 

e — ■ 

In the process of integrating over l x , two non-analytic 
contributions arise from (|D27|) . One comes from l x > 
(7^) 1 / 3 and goes like ^^1/3 , Plugging this back into 
and subtracting (and neglecting) a constant term, 
we obtain 



/+00 
-00 



-t-oo pQ 

dil / dui 

00 Jo 



\ix\ (iS(Q — u>) + iS(w) — vfIxY 



(D28) 



where we subtracted ±1,4(0,0) (hence the notation ^11^) 
and expanded in q. 

We further simplify the integrals, noticing that the 
fermionic propagator is dominated by vpl x since l x ~ 
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l/3\ 



, Nmg 2 



, Nmt; 2 
Wf7 5/3 Jh 



ft 



^273 
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where we substituted ~ ^fi 1 / 3 ) in the last steps. 

The other non-analytic contribution from 1|D27|) comes 
from typical v F l x ~ £(f2). It can be seen from an expan- 
sion of (|D27(1 for small values of both 8 and e, and goes 
like: 

(v F l x ) 2 + c 2 S(r!) 2 



2v 2 Fl \n\ 



■log((^y 2 + c 2 £(ft) 2 ).) (D30) 



One can explicitly verify that to get the logarithm, we 
only need the polarization operator 11(7, f2) to order 1/Z 3 . 
Like we said in Appendix C, to this order, polarization 
bubble can be evaluated with the full fermionic Green's 
functions but without vertex corrections. 

Substituting (|D30|I into the expression for 11^ and sub- 
tracting a constant part, we obtain: 



. 8Nmg 2 



2n 



(W 



(2ir) 5 jv F ' 1 J Q 
" cos Uv F l x ) 



+ OC 



00 dn 



dir. I 



(I 



(iE(fi — oj) + iTi(u) — v F l x ) 

x io g (M,) 2 + c 2 s(r!) 2 ) 

— lu) + iT,(lu) + v F q cos 9 — v F l x 

(D31) 

Using the 9 < — ► —9 symmetry, and splitting the inte- 
gral over O into two parts, one can rearrange this expres- 
sion as: 



n {2) ( q ,o) = i 



. Nmg 2 



—5 — 2I 



d6 



cos 9(1 + z 2 ) 



■ S(Q-^)+SM _ 

1 cS(n) 



log (l + 



; s(n-^)+s(Lj) 

cS(f2) 



1>F<? COS # 

cs(n) 
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where we defined z = v F l x /(cE(Q)). 

Let's now isolate the integral over z, given by: 



J = 



+°° (l + z 2 )log(l + z 2 ) 
-00 (ia — z) 3 (ia + 6 — .z) ' 
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where a = E(n ^+ S H and b = Sgg^ , and a > 0. 

Performing the contour integration in the lower half- 
plane (where lies the branch-cut), one gets: 

j--*, j 



-2tt 



dy 
dy 



(y + a) 3 (y + a-ib) 

{l-y 2 )(y + a + ib) 
(y + a) 3 ((y + a) 2 +b 2 )' 
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Once (|D34|) is plugged back into (|D32|) , only the imag- 
inary term survives due to the symmetry of the integral 
in 9. We are left with: 



n^M) 



ANmg 2 



2 /■"'/2 r+ca r+00 

q 2 d9 dy dVL 



C 2 -K^V F J ., 1 



dw- 



cos 2 9 



S (^) 2 [c^ ((I - w) 2 / 3 + W 2 / 3 ) + y] 

[y + c-i ((1 - ™) 2/3 + w 2 / 3 )] 2 + (2^) 

(D35) 



where we changed variables, defining w = lu/Q. 



Introducing the new variable t = ( cS ( fi ) ^ ^tris 

\ vfQ cos w / 3 



rewrites: 



ni 2) (s,o) 



4A r ^g 2 3/2 



0/2 4 3/2 l/2 y / 
(f/^TT^JVp LUq JO 



tt/2 



^q 3 / 2 / d9 (COS0) 3 / 2 



1 



(■+00 <>1 

x / dy dw- 

Ji Jo [c- 1 ({1 - w) 2 / 3 + w 2 / 3 ) + y] 

,+co j _ 2 

x / dt- 







1 + t 4 / 3 [ W + c- 1 ((1 - w) 2 / 3 + w 2 / 3 ] 



1 2 ' 
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A final change in variables leads to: 



n ^' )- 3/ 2 T'3/2 I7^ /2 / ^(-^) 3/2 
C 3 / Z TT^^V F W JO 

/" +0 ° , , !-y 2 

x / cty / 

Ji Jo [c- 1 ((1 - w) 2 / 3 + w 2 / 3 ) + y 



4Nmg 2 3/2 r-/ 2 



+00 
x / 



1 



1 + v^ 3 ' 



0/2 
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where v = t [y + cr 1 ((1 - w) 2/3 + w 2 / 3 )] 3/2 

Performing the integral over y, one is left with three 
independent integrals contributing to the numerical pref- 
actor: 



n^M) = — 3 ftj 2 1/2 ^ 2 r«> (™°) 3/2 

Z 1 5c+2((l- W ) 2 / 3 + W 2 / 3 ) 
x / dw ^j- 

(c+(l-w) 2 / 3 + w 2 / 3 ) 1 
f +00 ^ 

x I dv -j-. 

o 1 + v 4 / 3 



(D38) 
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These integrals can be performed separately and read: 



tt/2 



d9 (cos0) 3/2 



+ OC 



dv 



„ l + W 4 /3 

5c + 2s(w) 

dw ^ 

it (c + s(u>)) ' 



6[r(|)] 

4 : 
0.9438 



0.8740 



3.3322 
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and we used the notation s(w) = (1 — w) 2 / 3 + w 2 / 3 . 

Collecting all integrals, and rearranging the prefactor, 
the final result for the contribution of the first two dia- 
grams then writes: 



4 2) M) 



-0.1053 yJk F q 



,3/2 
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4. Other two diagrams 



The computation of the other two diagrams in Fig. 
proceeds along the same way. Far away from criticality, 
when 7i>f£ 2 /(1 + A) is small, and one can just expand 
perturbatively in the interaction, the sum of these two 
"drag" diagrams, which we label here and in the main 
text as Hb, was shown in [37j to be equal to 11a (g, 0) to 
the leading order in g (which in our model is g 3 , see (|D20I 
- ID22j) . Near criticality such simple relation no longer 
holds, but Ua (q, 0) and Hb(q, 0) remain of the same sign 
and of comparable magnitude. 

At criticality, we obtained for Ub(q> 0) 



n B («z,o) 



V2 



43 3 / 4 7T 4 
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where in rescaled variables (e.g., momentum is in units 
of q) 



1=1 dftdxx 



2tt 



d0- 



s 2 (x,n, 



, S 1 (x,Sl)S 1 (x + cos0,CL) 
and S(x, ft, 9) and Si{x, ft) are given by 



(D42) 



S{x,fi,9) = / doj 

Jo Jo 



de-L 



zS* — x cos(#i) 



1 



x cos t>i 



Si(x,(l) = / 
Jo 



(E* 
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where we introduced E* = E*(w,fi) = {fl- uj) 2 ^ 3 + lu 2 / 3 . 

We could not evaluate this integral explicitly, and we 
compute it under two simplifying assumptions 

• we compute Si(x, Q) by expanding to leading order 
in (E*/x) 2 , evaluating the frequency integral and 
plugging the result back into denominator. This 
way, we approximated Si(x,Q) by 



Si(x,n) 



n 



x 2 + (cVl 2 / 3 ) 2 



(D45) 



where c as 1.2 (see (|A6|) This procedure is similar 
to the one which led to (|A8I) . but here we cannot 
justify that only 1.x and 1/x 3 terms are relevant. 

• We replace E* by the same eft 2 / 3 in the integrand 
for S(x, Q, 9) 

• We assume that internal momenta are larger than 
the external one, i.e a; > 1 (x is measured in units 
of q), neglected terms 091) compared to 0{x) and 
set the lower limit of the integration over x at some 
number b. 

• We choose b by applying the same approximate 
computation scheme to IT4(q,u;) and requesting 
that the result coincide with the exact expression, 
Eq. Ip40| . 

Carrying out this calculation for II j B(g,0) we obtain 

n { 2\q,0) « -0.14 v/^f <Z 3/2 - (D46) 
This is the result that we cited in the text. 



APPENDIX E: TWO-LOOP RENORMALIZATION 
OF THE CHARGE SUSCEPTIBILITY 



In this Appendix we show that the singular contribu- 
tions to the static charge susceptibility from individual 
diagrams cancel out in the full expression of Xc(o)- The 
cancellation of the singularities in the charge responcc 
has been extensively studied in ID systems |41| . 



There are 10 different two-loop diagrams for the charge 
susceptibility, presented in Fig. [HJ The last two diagrams 



are identical to the ones we considered in the main text. 
We already argued there that these two diagrams cancel 



3G 



out in the case of a QCP in the charge channel. 

The other eight diagrams have to be considered to- 
gether. We demonstrate that the total contribution from 
these eight diagrams vanishes once one linearizes the dis- 
persion of the intermediate fermions. This still leaves out 
the contributions from non-linear terms in the dispersion, 
but one can show that these contributions are regular. 

To begin, consider one of these diagrams, e.g., diagram 
7 in Fig. [SJ In analytic form, the contribution from this 
diagram is: 

n 7 (g) = 2 J d 2 q 1 dui 1 d 2 q 2 duj 2 G k G k+q G k+qi G k+q2+q 

x G k + qi + q2 G k+qi+q2+q x qi Xq2i (El) 

where we labeled qi = (qi, LUi), and the combinatoric fac- 
tor 2 comes from the summation over spin indices. 
Introduce now 



G k G k 



1 



{G k - Gk+ qi ) , 



where 



qi cos ( 



(E2) 



(E3) 



and Qi is the angle between k « kp and q^. Shortening 
the notations further as qx = 1 and q 2 = 2, using the 
symbolic notation Jl „ for the 6-dimensional integral over 
qi and q 2 , and applying ljE2fl . we obtain: 



n 7 (<?) 

Similarly, 

n 8 (<z) = 2 

and further 
IL 5 (q) = 2 

n 3 (<?) 



2 / X1X2 

'1.2 



Gfc+qGfc+l-|-2 

2 2 

G k+q G k+ i 



X1X2 



XlX2 



1.2 



2 / X1X2 
'1,2 



ai(af - a 2 2 ) 

G k + q G k +i+ 2 
a\a 2 {a\ + a 2 ) 2 

G k+q G k+ i +2 ^ 
aia 2 (ai + a 2 ) 2 



Gfc+qGfe+l-|-2 

aia 2 (ai + a 2 ) 2 
(E4) 



, G k + q G k+ i 
'a\{a\ - a 2 ) 

(E5) 

J a\{a\ -a\) 

(E6) 



„2„,2 



+2- 



„2„,2 



n 2 (g) = 4 / X1X2- 

'1,2 



Gfc+qGfc+i 
2 2 

afo;| 



/a? 




W 








"1 





(E7) 



(E8) 



Mi) = * [ -2 



X 1X2 

2 *1 



Gfc+oGfc+1+2^ — " — o 

af - a| 
3a| -a? " 

-•-rfe+nLffc+l 2 5~ 

af — at 



(E9) 



114(g) — / ^^2 [Cfe+9Cfc+i+2 — 2Gfc +g Gfe + i] (E10) 

J 1,2 a l a 2 

and finally, 

n 6 (g) 



X1X2 

1,2 "l"! 



Gfe+qGfe+1^ 



«1 



(«i + a 2 ) 2 
— 2Gfc +? Gfe + i] (EH) 



Collecting the prefactors for Gfc+ g Gfc+i+ 2 and G k + q G k +i 
from all of the eight contributions we find that they can- 
cel out. 



APPENDIX F: MASS-SHELL SINGULARITY 

In this Appendix, we analyze in more detail the form of 
the self-energy near the fermionic mass shell. The inter- 
est to the mass-shell behavior of the self-energy was trig- 
gered by recent studies of the self-energy near a mass shell 
in a 2D Fermi liquid and for 2D Dirac fermions 
In both cases, the lowest-order self-energy diverges at the 
mass-shell, which forces to re-sum the perturbative series. 

At first glance, the same situation holds in our analy- 
sis at the QCP. Evaluating the self-energy in a two-loop 
expansion around free fermions and using the fermionic 
dispersion with the curvature, we obtain near the mass 
shell 



Aloe 



a- 



(Fl) 



This result implies that the "effective" quasiparticle 
residue for the Eliashberg theory, Z e ff cx dYj/de k log- 
arithmically diverges on the mass shell of free fermions. 
Without the curvature of the dispersion, the divergence 
would be stronger than logarithm. 

The issue we now have to address is whether Z still di- 
verges on the mass shell if we expand around the Eliash- 
berg solution, i.e., around fermions with 



G (*,w) 



1 



iS(w) - e k 



(F2) 



where, we remind, = lu + S(u;). 

It turns out that this is not the case: the expansion 
around the Eliashberg solution leads to a finite residue 
Z. At the two-loop order, we obtain, instead of IjFlf) 



log ; 

where 

ipz,x> = (1 



\2/3 



(1 



£fe) 



-') 2/3 + (z + z'-l 



,2/3 



(F3) 



(F4) 
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For simplicity, we restricted ourselves to the quantum 
critical regime where S(w) rs If ip g z i were equal to 

a constant, as it is when the system is in the Fermi liquid 
regime, and = Xu>, Z would diverge at u) = e/ (1+A). 
However, since S(f2 — w)+S(w) does not reduce to S(O), 
we have two additional integrations over z and z', and the 
logarithmic singularity is washed out. In particular, at 
€k = i.e. at the "Matsubara mass shell", we have 



in which case Z e ft is just a constant. Combining this 
with our earlier result that the renormalization of e& is 
also finite, Eq. H5.27|l . we obtain for the full fermionic 
Green's function at the smallest lu and 



G{k,w) 



Z, 



eff 



i£(w) - e* k 



(F6) 



J eff 



1 

N2 



Tog^7V-4.081og7V + 2. 



(F5) where e* k differs from ek by a constant factor. 
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